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I.  RESEARCH  OBJECTIVES 


The  overall  objective  of  the  research  is  to  develop  and  verify 
mathematical  models  of  aelamination  and  transverse  fracture  which  account  for 
local  (crack  tip)  and  global  damage  distributions.  One  specific  objective  is 
to  demonstrate  theoretical ly  and  experimentally  that  "work  potentials"  (which 
are  analogous  to  strain  energy)  exist  for  composites  with  constant  and 
changing  damage  and  with  viscoelastic  behavior;  validity  of  the  very  useful  J 
integral  theory  of  fracture  analysis  depends  on  the  existence  of  work 
potentials.  The  second  specific  objective  is  to  develop  and  verify  methods  of 
analysis  for  predicting  crack  growth  in  elastic  and  viscoelastic  composites 
with  distributed  damage  using,  when  applicable,  the  J  integral. 

II.  RESEARCH  SUMMARY 


1.  Overview 

Methods  of  deformation  and  fracture  characterization  and  prediction  are 
simplified  when  strain  energy-like  potentials  based  on  mechanical  work  can  be 
used,  ^s  described  in  the  first  paper  in  the  Appendix,  "Deformation  and 
Fracture  Characterization  of  Inelastic  Composite  Materials  Using 
Potentials".  With  these  so-called  work  potentials,  important  theoretical  and 
experimental  methods  using  the  J  integral  and  energy  release  rate  (originally 
developed  for  fracture  of  elastic  media  and  fracture  Initiation  in  metals  with 
plastic  deformations)  may  be  extended  to  fracture  initiation  and  crack 
propagation  in  monolithic  and  composite  materials. 

Section  2  is  concerned  primarily  with  a  study  of  rectangular  angle-ply 
bars  under  combined  axial  and  torsional  deformation.  This  investigation  is  by 
Mark  Lamborn,  a  Ph.D.  student,  and  it  deals  with  (i)  the  existence  and  then 
determination  of  a  work  potential  for  specimens  with  significant  amounts  of 
distributed  damage  and  (1i)  characterization  and  prediction  of  delamination 
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when  distributed  damage  and  the  mode  III  (antiplane)  component  of  energy 
release  rate  are  relatively  large.  Experimental  studies  to-date  on  a 
composite  with  brittle  resin  (Hercules  AS4/3502),  as  reported  in  the  first 
paper  in  the  Appendix,  and  on  a  composite  with  rubber-toughened  resin  (Hexcel 
T2CT145/F155) ,  as  reported  in  Section  2,  show  that  a  work  potential  exists 
with  proportional  axial  and  torsional  deformation.  On  the  basis  of  the 
theory,  we  expect  a  work  potential  to  exist  over  many  sets  of  non-proportional 
deformation  histories;  some  experimental  work  planned  for  the  near  term  on  a 
new  AFOSR  grant  will  be  concerned  with  non-proportional  deformation. 

Also  sumnarized  in  Section  2  is  work  on  cyclic  axial-torsional  loading  of 
tubes  by  another  Ph.D.  student,  Richard  Tonda.  Here,  the  objective  is  to 
develop  a  method  for  removing  from  experimental  data  effects  of  visco¬ 
elasticity  on  hysteresis  and  stiffness  so  that  effects  of  distributed  damage 
growth  may  be  more  readily  studied. 

A  numerical  investigation  of  crack  growth  initiation  and  continuing  crack 
growth  in  inelastic  media  is  sunmarized  in  Section  3;  this  work  was  done  by 
Randall  Weatherby,  a  Ph.D.  student.  A  multivalued  work  potential  is  used  to 
characterize  material  behavior.  Path  independence  of  the  J  integral  is  then 
studied  and  shown  to  be  consistent  with  theory  developed  on  an  earlier  AFOSR 
grant. 

Finally,  in  Section  4  a  delamination  study  by  Douglas  Goetz,  a  Ph.D. 
student,  is  described.  Mode  I  delamination  of  laminates  with  multiple  fiber 
orientations  is  investigated,  and  the  usefulness  of  the  J  integral  as  a 
fracture  characterizing  parameter  is  demonstrated. 

Most  of  the  work  described  in  Sections  2  and  4  has  not  yet  been 
published,  and  therefore  detailed  summaries  are  given. 
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2.  Studies  of  Laminates  Under  Axial  and  Torsional  Deformation 

The  research  effort  last  year  was  directed  toward  selection  of 
rectangular  bar  specimens,  development  of  experimental  procedures,  and 
development  of  data  reduction  schemes.  The  objective  in  specimen  selection 
was  to  determine  which  specimens  which  would  exhibit  a  significant  coupling 
effect  of  rotation  and  axial  deformation  on  loads  in  order  to  be  able  to 
critically  evaluate  the  work-potential  theory.  Laminates  of  various  layups 
and  "'ometrles  were  tested.  Test  specimens  of  Hercules  AS4/3502 
graphite/epoxy  were  24  plies  thick  and  either  balanced  angle-ply  or  other 
balanced  symmetric  laminates.  Tests  were  run  under  conditions  of  proportional 
and  non-proportional  straining;  instrumentation  difficulties  voided  the 
results  for  non-proportional  straining  and,  in  fact,  resulted  in  a  delay  of 
testing  for  several  months.  All  layups  and  geometries  were  found  to  exhibit 
only  a  small  effect  of  rotation  on  axial  load,  while  the  effect  of  axial 
displacement  on  torque  was  significant.  Test  data  from  proportional  straining 
tests  (i.e.  the  twist  was  proportional  to  the  axial  extension)  were  evaluated 
for  the  existence  of  a  work  potential  as  discussed  in  [2.1];  references  for 
section  2  start  on  page  22.  The  results  of  some  of  these  tests  appear  In 
Figures  3  and  4  of  [2.1],  which  is  in  the  Appendix.  These  results  provided  at 
least  a  limited  check  of  the  existence  of  a  work  potential  for  conditions  of 
proportional  straining. 

Experimental  work  was  continued  in  the  current  reporting  year  to  verify 
the  existence  of  a  work  potential  for  fiber-reinforced  plastic  laminates 
subjected  to  combined  axial  and  torsional  loadings.  Hexcel  T2CT145/F155,  a 
material  system  consisting  of  a  rubber  toughened  epoxy  matrix  with  graphite 
fibers,  was  used  to  fabricate  [ ±35 ° ] gs  flat  bar  specimens.  These  were  tested 
under  the  conditions  of  proportional  straining.  Typical  test  results 
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expressed  in  terms  of  "nominal"  stresses  and  strains  appear  as  the  dotted 
lines  in  Figures  2.1  and  2.2.  The  nominal  stresses  and  strains  are  defined  by 
Equation  (39)  in  (2.1).  The  test  data  was  evaluated  for  the  existence  of  a 
work  potential  by  the  method  outlined  in  [2.1].  The  "nominal"  axial  stress- 
strain  responses  for  cases  of  combined  axial  and  torsional  loadings,  were 
predicted  from  the  experimentally  determined  work  potential.  The  predicted 
responses  are  shown  as  the  solid  lines  appearing  in  Figure  2.2.  The 
predictions  shown  in  this  Figure  represent  tests  run  at  intermediate  values 
of  e/y,  the  ratio  of  nominal  axial  strain  to  shear  strain.  These  tests  span 
the  levels  of  e  and  y  where  sufficient  data  were  available  for  calculation 
purposes.  Tests  were  run  at  c/y  =  .25,  2,  and  4  but  could  not  be  predicted 
since  they  represent  limit  cases  and  a  lack  of  data  exists  at  the  strain 
levels  reached  in  these  tests.  As  seen  in  Figure  2.2  the  predicted  responses 
agree  well  with  experiment.  It  is  tentatively  concluded  that  a  work  potential 
exists  for  conditions  of  proportional  straining  and  increasing  material 
damage. 

Data  from  the  proportional  straining  tests  on  laminates  constructed  of 
the  Hexcel  material  system  were  analyzed  to  determine  an  analytical 
representation  of  the  work  potential.  The  theory  in  [2.1]  was  used  to 
characterize  the  material  and  make  predictions  of  material  response  for 
conditions  of  combined  axial  and  torsional  loadings.  The  relationship  used  to 
characterize  the  test  data  is  that  corresponding  to  Equation  (28)  of  [2.1], 
The  form  used  in  this  work  is 


WT  =  W  +  WQ  (2.1) 

where  Wy  =  the  total  mechanical  work,  W  =  the  strain  energy  density  function, 
and  Wq  =  the  work  of  damage.  All  quantities  appearing  in  Equation  (2.1)  are 
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expressed  on  a  per  unit  initial  volume  basis. 

The  total  mechanical  work  input  is  defined  by 


“t  -  { 


(2.2) 


where  qi  =  the  generalized  displacements  and  Qi  =  the  generalized  forces 
(i=l,2).  The  integral  is  taken  from  generalized  displacement  state  1  to 
generalized  displacement  state  2.  For  stable  microcracking  and  some  other 
damaging  mechanisms,  it  has  been  shown  that  the  integral  appearing  in  Equation 
(2.2)  is  independent  of  path  [2.1].  Wj  will  be  referred  to  as  the  work 
potential .  From  Equation  (2.2) 


> 

Qj 

=  aWy/sqj 

(2.3) 

i 

A  basic 

assumption  of 

this  theory  is  that 

the 

strain  energy  density 

j.' 

function 

W  -  W(q1fDk) 

exists  and  is  defined  by 

s 

qj 

=  sW/sq^. 

(2.4) 

1] 

where  Dk 

(k»l  ,2 . K) 

are  damage  parameters. 

The 

set  Dk  includes  all  the 

parameters  needed  to  account  for  any  structural  changes  occurring  in  the 
material  during  loading.  Such  structural  changes  could  include  microcracking, 
delamination,  and  fiber  breakage.  For  the  purposes  of  obtaining  an  initial 
analytical  representation  of  the  experimental  data  it  is  assumed  that  one 
damage  parameter  accounts  for  all  the  structural  changes.  This  damage 
parameter  is  taken  to  be  Wg.  The  damage  Wg  will  change  with  the  applied 
displacements  q^;  thus  Wg  =  Wg(qi).  From  Equations  (2.1),  (2.3),  and  (2.4) 
the  equation  for  predicting  the  current  value  of  Wg  becomes 


3W/3Wg  =  -1 


(2.5) 


* 


8 


provided  this  equation  predicts  dWg/dt  >  0.  For  cases  when  dWg/dt  <  0  it  can 
be  shown  from  thermodynamics  that  Wg  must  be  constant. 

For  the  test  data  under  examination  the  generalized  forces  and 
displacements  are  taken  as  q^  =  e,  q^  =  y.  Qj  =  o,  =  t/3  where  e,  y.  o.  and 

t  are  the  nominal  strains  and  stresses  defined  in  Equation  (39)  of  [2.1].  In 
terms  of  these  variables  Equation  (2.2)  is  rewritten  as 


WT  =  f  a  de'  +  J  i  dr' 
0  0  J 


(2.6) 


where  the  initial  values  of  the  generalized  displacements  are  taken  as  zero 
(unstrained  state).  For  the  purposes  of  analyzing  the  experimental  data  the 
strain  energy  density  function  is  taken  to  be 


u  1  r  2  1  r  2 
W  -  75-  C  e  +75C  y 

CO  C  I 


(2.7) 


where  Cq  and  are  functions  of  Wg.  The  strain  energy  density  function,  W, 
as  given  by  Equation  (2.7)  is  that  of  a  linear  elastic  material  when  Wg  is 
constant.  Substitution  of  Equations  (2.6)  and  (2.7)  into  Equation  (2.1) 


yields 


wn  N  0  dt'  +  J  3 
u  0  0  J 


f  1  dr'  -i  C 


1  f  2 
2  Cx  Y 


(2.8) 


from  which  Wg  is  seen  to  be  the  sum  of  the  shaded  areas  under  the 
a  versus  c  and  t/3  versus  y  curves  as  shown  in  Figures  2.3  and  2.4.  From 
Equation  (2.4) 


Area  under  the  axial  stress-strain  curve  contributing  to  the  work  of 
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Equation  (2.8)  shows  that  when  W  is  evaluated  on  the  loading  curve,  Cq  and 
correspond  to  the  secant  moduli  of  the  a  versus  e  and  x/3  versus  y  curves, 
respectively. 

Equations  (2.8)  and  (2.9)  were  used  to  obtain  values  of  Cq  and  as 
functions  of  Wp  for  proportional  straining  tests  run  at  e/y  =  0,  .5,  .75,  1, 
1.5,  2,  4,  and  °°.  The  results  of  these  calculations  appear  in  Figure  2.5.  It 
was  found  that  the  test  ran  at  e/y  =  1.5  experienced  the  greatest  amount  of 
damage  as  defined  by  the  maximum  value  of  Wq.  Analytical  expressions  were 
then  obtained  for  C  and  C  as  functions  of  Wn  for  this  one  test.  These  are 
shown  as  the  solid  lines  in  Figure  2.5.  The  analytical  expressions  for  these 
curves  were  used  with  Equations  (2.5)  and  (2.9)  to  predict  the  material 
response  for  tests  run  at  the  same  values  of  e/y  as  those  appearing  in  Figures 
2.1  and  2.2.  The  results  of  these  predictions  are  shown  in  Figures  2.6  and 
2.7.  These  figures  show  that  the  predicted  responses  for  e/y  >  1  agree 
reasonably  well  with  experiment  while  those  for  e/y  <  1  do  not.  The  reason 
for  this  appears  to  be  that  a  second  damage  mechanism  becomes  significant  at 
the  larger  levels  of  twist.  This  would  require  that  a  second  damage  parameter 
be  incorporated  into  the  model.  The  additional  damage  mechanism  is  apparent 
for  the  data  shown  in  Figure  2.8.  Appearing  in  Figure  2.8  are  the  secant 
moduli  plotted  against  each  other  for  e/y  =  .5,  .75,  1,  1.5,  2,  and  4.  These 
data  show  that  the  curves  for  the  different  tests  are  essentially  the  same 
except  for  e/y  =  .5  and  .75;  the  two  tests  correspond  to  the  greatest  levels 
of  twist  reached  for  conditions  of  combined  axial  and  torsional  loadings. 
Characterization  of  the  second  damage  mechanism  is  presently  under  study. 

Physically,  the  second  damage  mechanism  may  be  delamination.  For  the 
angle-ply  layup,  delamination  would  have  a  large  affect  on  the  torsional 
stiffness  but  only  a  small  affect  on  the  axial  stiffness. 


-ted  and  experimental  shear  stress-strain  curves  with  proportional  axial  and  tor- 
1  straininq.  Hexcel  Fl 55;  24  plies  [±35].,;  8/5  in.  long  x  0.5  in.  wide  x  0.145  in 


Predicted  and  experimental  axial  stress-strain  curves  with  proporti 
torsional  straining  for  the  same  tests  used  for  Figure  2.5. 


Curves  of  CT  versus  C  showing  the  change  in  material  behavior  at  large  torque 
levels  due  to  second  damage  mechanism. 


New  Test  Methods  and  Equipment:  During  the  course  of  these  tests, 
experimental  procedures  were  refined.  Deficiencies  in  measurement  procedures 
were  recognized  and  corrected.  A  computer  program  was  developed  to  check  data 
from  a  series  of  proportional  straining  tests  for  the  existence  of  a  work 
potential.  The  program  was  tested  using  data  from  a  series  of  tests  on 
aluminum  under  combined  axial  and  torsional  loads.  The  program  verified  the 
existence  of  a  work  potential  for  the  aluminum  well  into  the  yielded  state,  as 
expected. 

New  equipment  has  been  purchased  to  aid  in  the  experimental  work.  A  load 
cell  which  is  more  sensitive  to  low  torque  levels  than  the  one  used  last  year 
has  been  acquired  and  is  currently  in  use.  The  fiber-reinforced  plastic 
laminates  that  have  been  tested  in  combined  axial  and  torsional  loading 
typically  experience  maximum  torque  levels  of  350  in-lbs.  The  new  load  cell 
has  torque  range  settings  of  400  in-lbs.  and  200  in-lbs,  while  the  previous 
load  cell  had  a  minimum  torque  range  setting  of  1000  in-lbs.  The  axial 
capacity  of  both  old  and  new  load  cells  is  the  same. 

A  new  computer  system  was  purchased  to  aid  in  data  acquisition  and 
reduction.  The  new  system  provides  increased  data  acquisition  speed  as  well 
as  increased  data  handling  capabilities.  This  system  is  currently  being  used 
in  the  experimental  work. 

Problems  were  experienced  with  the  computer  controller  for  the  MTS  axial- 
torsion  testing  machine.  The  problem  involved  the  unplanned  loading  of  test 
specimens  when  control  switched  from  manual  to  the  computer.  In  some 
instances  this  damaged  or  failed  the  test  specimens.  Computer  control  is 
beneficial  when  complicated  displacement  and  rotation  histories  are  used  as 
input.  Work  is  currently  underway  by  MTS  to  resolve  the  problem. 

Related  Studies:  A  literature  survey  was  performed  to  review  prior  work 


v?  ■r.V.'vv'-- 


concerned  with  mode  III  fracture  of  metals  and  composite  materials.  It  was 
found  that  little  data  were  available  on  mode  III  fracture  characterization  of 
materials.  The  majority  of  the  available  data  were  from  studies  on  metals. 
An  investigation  of  the  experimental  methods  used  to  characterize  metals 
subjected  to  mode  III  deformation  indicates  that  these  methods  would  be 
unsuitable  for  mode  III  delamination  studies  of  f iber-reinforced  plastic 
laminates. 

Studies  concerned  with  the  mode  III  fracture  characterization  of  fiber- 
reinforced  composite  materials  are  few.  Experimentally  determined  values  for 
the  critical  energy  release  rate,  Gjjjc,  have  been  obtained  for  a  randomly 
oriented  short  fiber  composite  [2. 2|  and  a  fiber-reinforced  plastic  laminate 
[2.31.  In  both  cases  the  measured  value  of  Gj j jc  was  found  to  be  greater  than 
measured  values  of  Gjc  and  Gjjc  for  the  same  materials.  For  the  case  of  the 
fiber-reinforced  plastic  laminate  the  critical  energy  release  rate  in  mode  III 
was  determined  for  three  different  layups.  The  value  of  G  j  j Ic  was  found  to  be 
greater  than  GIc  and  G  j  j  c  for  all  layups.  A  recent  finite  element  study 
determined  energy  release  rates  for  laminates  subjected  to  pure  torsional 
loading  [2.4].  For  the  three  layups  studied,  the  energy  release  rate  in  mode 
III  was  found  to  be  greater  than  that  in  modes  I  and  II,  and  pure  mode  III  was 
achieved  in  a  [0°lgs  laminate. 

Analytical  Studies  Related  to  Torsion  and  Torsion-Induced  Delamlnatlon: 
To  gain  a  better  understanding  of  the  torsion  problem  of  a  laminate,  a  closed 
form  solution  to  the  linear  elastic  torsion  problem  for  a  cross-ply  laminate 
was  obtained.  The  solution  satisfies  the  field  equations  on  a  ply-by-ply 


basis.  Displacement  and  stress  continuity  conditions  are  required  at  the 
interfaces  between  0°  and  90°  plies.  The  result  is  a  Fourier  series  solution 
for  the  stress  function  of  each  ply.  Convergence  of  the  solution  at  the 


interfaces  between  plies  is  not  necessarily  guaranteed  and  is  currently  under 
investigation.  This  type  of  solution  may  provide  reasonable  predictions  of 
material  response  for  angle-ply  laminates  with  a  large  number  of  plies  and  may 
be  used  in  test  specimen  design. 

The  torsion  problem  of  cracked  isotropic  plates  was  also  studied 
theoretically  and  experimentally.  An  isotropic  plate  with  a  through -the- 
thickness  crack  was  analyzed  (Figures  2.9  and  2.10).  Shown  in  Figure  2.9  is 
such  a  plate  subjected  to  torsional  loading.  As  indicated  in  Figure  2.9,  the 
separated  sides  of  the  plate  will  bend  as  beams  under  torsional  loads.  It  was 
an  important  aspect  of  this  study  to  examine  the  effect  of  grip  restraint  on 
bending  and  thereby  determine  the  increase  in  torsional  rigidity  of  the 
separated  plate.  The  torsional  rigidity  in  the  presence  of  the  crack  (Figure 
2.10)  was  predicted  using  standard  torsion  theory  and  a  strength  of  materials 
approach  which  accounts  for  the  bending  of  the  separated  sides.  Figure  2.10 
shows  the  results  of  a  test  made  on  an  aluminum  plate  and  the  predictions  made 
by  the  two  theories.  It  is  apparent  from  these  results  that  resistance  to 
bending  has  a  significant  affect  on  the  torsional  response  of  separated 
plates.  These  results  suggest  that  bending  resistance  of  the  grips  may  also 
have  a  significant  effect  on  the  torsional  response  of  fiber-reinforced 
plastic  laminates  with  edge  delaminations. 

Additional  Studies  of  Rectangular  Bars:  Research  on  a  new  grant  will 
concentrate  on  two  areas  needing  further  investigation.  One  topic  area  deals 
with  the  question  of  the  existence  of  a  work  potential.  The  experimental  work 
to-date  has  verified  the  existence  of  a  work  potential  under  conditions  of 
proportional  deformation  when  the  deformations  increase  in  time.  The 
existence  of  a  potential  for  these  conditions  does  not  guarantee  its  existence 
for  all  deformation  histories.  Future  work  will  investigate  the  existence  of 


20 


experimental  and  predicted  torque  reponse  curves. 


a  work  potential  for  non-proportional  deformation  histories,  including 
histories  where  the  deformations  decrease  in  time.  These  input  deformation 
histories  should  include  processes  where  some  damage  mechanisms  would  remain 
inactive  while  others  would  be  active.  These  would  correspond  to  different 
branches  of  the  multivalued  potential.  The  existence  of  a  potential  along 
these  different  branches  needs  verification.  The  question  of  the  existence  of 
a  work  potential  when  viscoelastic  effects  are  significant  also  needs 
investigation.  For  such  cases  as  these,  the  theory  outlined  in  [2.1]  and 
(2.5[  will  be  used  to  check  for  the  existence  of  a  work  potential.  For 
processes  where  it  is  verified  that  a  work  potential  exists,  the  generalized  J 
integral  theory  given  in  [ 2 . 5 J  will  be  used  in  mode  III  and  mixed  mode 
fracture  studies. 

Mode  III  fracture  characterization  of  fiber-reinforced  plastic  laminate 
is  an  area  requiring  further  investigation.  As  previously  mentioned,  studies 
of  this  type  are  few.  The  tentative  test  specimen  which  will  be  used  for  mode 
III  characterization  studies  is  a  rectangular  plate  laminate  with  pre¬ 
fabricated  edge  and  through-width  delaminations.  It  is  the  intent  of  future 
research  to  use  this  test  specimen  to  study  mixed-mode  and  mode  III 
delamination  of  fiber-reinforced  plastics  laminates.  At  present  a  [O0^ 
plate  of  the  Hexcel  composite  has  been  fabricated  with  teflon  release  paper 
used  to  form  pre-existing  edge  delaminations. 

A  X-ray  machine  is  being  purchased  so  that  X-ray  fractography  may  be  used 
to  monitor  delamination  and  other  damage  growth.  This  equipment  should  also 
be  useful  in  identifying  other  damage  mechanisms. 

Studies  of  Thin-Walled  Tubes:  A  Ph.D.  dissertation  involving  experimental 
work  on  tubular  specimens  was  completed  by  Richard  Tonda  during  the  grant 
period.  A  significant  part  of  the  study  was  concerned  with  development  of  a 
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technique  for  modifying  stress-strain  data  on  fibrous  composites  so  that  the 
effects  of  changing  damage  may  be  observed  without  the  complicating  effects  of 
matrix  viscoelasticity.  The  method,  which  is  based  on  micromechanical 
considerations  and  is  not  limited  to  tubular  specimens,  reduces  the  behavior 
to  that  of  an  equivalent  elastic  composite  with  damage.  The  fibers  are 
assumed  to  be  continuous  and  linearly  elastic.  The  theoretical  basis  is 
developed  and  then  the  method  is  illustrated  using  results  from  cyclic  axial- 
torsional  loading  of  tubular  specimens  of  graphite/epoxy  laminates.  This  work 
is  reported  in  12.6] . 
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A  Ph.D.  dissertation  on  analysis  of  crack  growth  was  completed  by  Randall 
Weatherby  during  the  grant  period.  The  work  is  summarized  here. 

In  most  materials,  macrocrack  extension  is  accompanied  by  inelastic 
phenomena  (such  as  microcracking  or  plastic  deformation)  throughout  a  region 
surrounding  the  crack  tip.  Immediately  ahead  of  that  crack  tip,  strain 
localization  occurs  in  a  small  volume  of  heavily  damaged  material  referred  to 
as  the  failure  zone  or  fracture  process  zone.  In  this  study,  the  failure  zone 
and  the  surrounding  zone  of  inelastic  material  are  treated  as  two  distinct 
regions.  The  failure  zone  is  assumed  to  be  thin  relative  to  its  length  and  is 
represented  in  a  two-dimensional  finite  element  model  as  tractions  which  act 
across  the  crack  faces  near  the  tip.  An  opening  mode  of  crack  tip  deformation 
is  assumed.  The  normal  traction  at  any  point  on  the  crack  surface  in  the 
failure  zone  is  specified  as  a  decreasing  function  of  the  crack  opening 
displacement  which  vanishes  after  a  critical  value  of  displacement  is 
reached.  Two  different  rate- Independent,  inelastic  continuum  characteri¬ 
zations  based  on  multivalued  work-potentials  are  used;  one  models  metal 
plasticity  and  another  represents  microcracking  in  brittle  materials.  Both 
constitutive  models  allow  for  the  definition  of  a  generalized  J  integral 
developed  by  Schapery,  which  has  the  same  value  for  most  paths  around  the 
crack  tip  for  realistic  distributions  of  plasticity  or  damage  in  the  material 
surrounding  a  stationary  or  propagating  crack.  This  path  independence  and  the 
equivalence  between  J  and  the  work  input  to  the  last  ligament  of  material  in 
the  failure  zone  are  verified  numerically  in  a  transient  crack  growth  problem; 
both  initiation  and  propagation  are  studied  under  conditions  of  small-scale 
inelasticity.  In  addition,  steady-state  crack  growth  is  studied  in  two 
different  specimen  geometries.  Simplified  J  integral  analyses  are  used  to 
estimate  the  work  input  to  the  failure  zone  for  these  steady-state  problems. 
The  J-integral  estimations  are  compared  with  finite  element  results  to 
determine  the  accuracy  of  the  simplifed  analyses. 


4.  Determination  of  the  Mode  1  Delamination  Toughness  of  Multidirectional 
Laminates 

The  objective  of  this  study  is  to  develop  a  procedure  for  determining  the 
mode  I  delamination  fracture  toughness  of  multidirectional  composite 
laminates.  To  date,  the  focus  of  the  work  has  been  on  application  of  the  J 
integral  to  double  cantilever  beam  (DCB)  tests  in  order  to  investigate  the 
importance  of  geometry,  layup,  rate  effects,  and  fracture  morphology.  The 
general  goal  is  to  synthesize  principles  which  can  guide  the  prediction  of 
delamination  performance  using  limited  material  data.  The  following 
discussion  highlights  the  method  used  and  results  obtained. 

The  material  used  in  this  study  is  T2CT145/F155,  manufactured  by 
Hexcel .  It  uses  a  rubber-toughened  epoxy  matrix  with  about  6  %  rubber  by 
weight.  The  nominal  fiber  volume  fraction  for  all  the  material  tested  was 
57%.  This  material  system  was  chosen  in  order  to  accentuate  possible 
nonlinear  and  viscoelastic  responses  for  a  commercially  available  toughened 
composite.  It  is  also  currently  being  used  for  other  work  on  this  contract 
and  other  contracts  at  Texas  A&M,  allowing  sharing  of  material  data.  Three 
layups  were  used.  Hereafter,  they  will  be  referred  to  using  the  following 


designations: 

Designation _ Layup _ No.  Plies 

unidirectional  (C^l  24 

fiber-dominated  [+-45/0g/-+45] anti Sym.  24 

angle-ply  [+-45/(-+45)2/(+-45)2/-+45]antisynK  24 

These  designations  will  also  be  used  in  the  captions  of  the  figures  to  be 

discussed  below.  The  unidirectional  layup  is  the  one  which  has  been  most 

commonly  used  in  the  literature  to  characterize  delamination.  In  this  study 

it  was  used  not  only  to  examine  delamination  between  zero  degree  plies  where 

ply  interpenetration  could  occur,  but  to  give  a  baseline  for  evaluating  data 
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for  multidirectional  layups.  The  fiber-dominated  layup  was  used  to  study 
delamination  at  a  +-45  degree  interface.  Each  leg  of  the  DCB  specimen  was 
made  balanced  and  symmetric  in  order  to  eliminate  stretching-shearing  and 
stretching-bending  coupling.  Antisymmetry  about  the  midplane  of  the  laminate 
was  designated  in  order  to  put  the  delamination  plane  at  a  +-45  degree 
interface.  The  stacking  sequence  for  the  angle-ply  was  chosen  to  minimize 
bend ing-twi sting  coupling  as  well. 

All  data  shown  in  the  following  figures  were  generated  using  DCB 
specimens.  The  fracture  toughness  was  calculated  using  the  J  integral  method 
of  [4.1],  which  allows  for  nonlinear  elastic  or  inel astic  behavior.  In  this 
formulation  the  J  integral  is  calculated  using  the  moment  at  the  crack  tip 
during  crack  extension  and  the  moment-curvature  relationship  for  one  leg  of 
the  specimen.  The  J  integral  is  twice  the  area  to  the  left  of  the  moment- 
curvature  curve  divided  by  the  specimen  width,  as  shown  in  Fig.  4.1.  In  the 
following  discussion,  J  will  correspond  to  this  integral  definition.  The 
energy  release  rate  was  also  calculated  for  comparison.  The  symbol  G  will 
signify  the  energy  release  rate  calculated  by  the  area  method  [4.2],  which 
allows  for  nonlinear  elastic  behavior.  The  area  method  comes  from  the 
derivative  definition  of  energy  release  rate,  and  provides  an  average  value 
for  an  increment  of  crack  growth.  The  symbol  G-j  will  refer  to  the  energy 
release  rate  calculated  using  the  equation  for  a  linear  elastic  beam. 
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where  Pc  is  the  load  at  crack  extension,  is  the  load  line  displacement  at 
crack  extension,  b  is  the  specimen  width,  and  a  is  the  crack  length. 

Motivation  for  using  a  J  integral  analysis  came  from  the  need  for  an 
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Figure  4.1  Double  cantilever  beam  specimen,  integration  path,  and  moment- 
curvature  diagram  for  calculation  of  J  [after  1] 


analysis  which  would  allow  for  continuum  damage  such  as  might  be  generated 
throughout  the  legs  of  a  DC8  specimen  during  testing;  the  value  of  0  is  only 
the  work  input  to  the  crack  tip  during  self-similar  crack  growth.  The  area 
method  of  calculating  the  energy  release  rate  allows  for  geometric  and 
material  nonlinearity,  as  does  J,  but  does  not  differentiate  between  energy 
which  goes  into  driving  the  delamination  and  that  which  causes  general  damage 
away  from  the  crack  tip,  leading  to  overestimation  of  toughness.  The 
difference  is  illustrated  in  Fig.  4.2,  where  G  and  J  are  plotted  versus  crack 
length  for  three  angle-ply  DCB  specimens.  Note  that  for  each  specimen  G  is 
significantly  higher  than  J.  The  presence  of  damage  (and  a  small  amount  of 
viscoelastic  behavior)  is  illustrated  by  the  moment-curvature  relationship  for 
the  angle-ply  layup,  as  shown  in  Fig.  4.3.  Note  that  the  unloading  portion  of 
the  moment-curvature  plot  does  not  retrace  the  loading  portion. 

For  layups  containing  a  high  percentage  of  0  degree  fibers,  one  would 
expect  that  G  and  J  would  give  better  agreement,  since  the  effect  of  continuum 
damage  on  beam  deformations  would  be  small.  The  moment-curvature 
relationships  for  the  unidirectional  and  fiber-dominated  layups  were  linear, 
and  the  unloading  portions  of  the  curves  retraced  the  loading  portions.  When 
this  is  the  case  and  linear  beam  theory  applies,  it  can  be  shown  that  J  is 
equal  to  G^  calculated  using  equation  1.  This  formulation  is  helpful  when 
comparing  energy  release  rate  and  0  because  G-j  is  calculated  for  a  particular 
crack  length  (as  J  is  for  the  linear  and  nonlinear  cases),  avoiding  averaging 
over  crack  lengths  as  in  the  area  method.  Figure  4.4  is  the  load-displacement 
record  from  a  fiber-dominated  specimen  which  displayed  the  development  and 
breakdown  of  a  tie  zone,  causing  the  delamination  resistance  to  increase  and 
then  suddenly  drop  with  crack  growth.  The  corresponding  G-j  and  J  are  plotted 
in  Fig.  4.5.  Both  initiation  and  arrest  values  are  shown.  The  two  methods 
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Figure  4.4  Load  vs.  displacement  plot  for  a 
The  layup  is  [+-45/0g/(-+45)2/08 


Figure  4.5  6^  and  J  vs.  crack  length  for  the  fiber-dominated  specimen  of 

figure  4.4.  Arrest  values  are  marked  by  squares  around  data 
points.  The  layup  is  |+~45/0ft/(-+45W0o/-+45| . 


are  in  reasonably  close  agreement,  supporting  the  interpretation  of  the 
difference  between  energy  release  rate  and  J  for  the  angle-ply  layup  as  being 
caused  by  damage.  Even  when  damage  is  not  present,  the  simple  expression  in 
equation  1  does  not  always  hold;  when  geometric  nonlinearity  is  present  due  to 
large  rotations  of  the  legs,  the  energy  release  rate  must  be  calculated  using 
either  the  area  method  (possibly  losing  information  due  to  averaging)  or  using 
nonlinear  beam  theory  (4.31. 

The  usefulness  of  J  as  a  characterizing  parameter  depends  on  its 
independence  of  geometry,  such  as  specimen  width.  To  investigate  the  effect 
of  specimen  width,  angle-ply  specimens  1/2,  1,  and  2  inches  wide  were 
tested.  Representative  results  are  given  in  Fig.  4.6.  The  high  values  for 
the  1  inch  wide  specimen  for  shorter  crack  lengths  were  found  to  be  associated 
with  the  complex  fracture  morphology  which  developed  near  the  starter  crack. 
For  more  nearly  self-similar  crack  advance,  it  appears  that  width  does  not 
have  much  effect  on  the  fracture  toughness  measured  over  the  range  of  widths 
and  crack  lengths  studied.  Recent  work  at  Texas  A&M  [4.4]  indicates  that  the 
state  of  stress  in  the  legs  of  the  DCS  specimen  is  in  the  transition  range 
between  plane  strain  and  plane  stress  over  the  practical  range  of  crack  length 
to  width  aspect  ratios.  However,  the  effects  of  the  transition  on  J  were  not 
seen  for  the  aspect  ratios  in  our  tests.  There  is  another  type  of  specimen 
behavior  which  does  show  a  marked  width  effect.  Unlike  unidirectional 
composites,  multidirectional  composites  exhibit  significantly  curved  crack 
fronts  in  DCB  tests.  The  curvature  depends  on  the  width  of  the  specimen. 
Figure  4.7  compares  the  normalized  crack  front  profile  for  1/2  inch,  1  inch, 
and  2  inch  wide  specimens  for  several  crack  lengths.  The  thickness  of  one  leg 
is  0.062  inch.  Apparently  the  crack  front  curvature  has  two  sources;  the 
anticlastic  curvature  of  the  DCB  leg  and  free  edge  effects.  The  general  trend 
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4.7  Normalized  crack  front  profiles  for  three  angle-ply  DCB 

specimens  of  different  widths.  The  profiles  for  various  crack 
lengths  have  been  superimposed. 
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Figure  4.7  (continued) 
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Figure  4.8  J  vs.  crack  length  for  three  different  layups 


(LB/IN) 


39 


measured.  Note  that  once  a  more  or  less  uniform  fracture  morphology  was 
established,  the  toughness  values  settled  into  better  agreement.  The  two 
angle-ply  specimens  shown  were  next  to  each  other  in  the  plate  of  material 
before  cutting.  Dark  and  light  bands  were  formed  when  the  crack  tip  jumped 
forward,  then  arrested.  Dark  regions  were  formed  during  rapid  crack  advance; 
light  regions  correspond  to  slow  crack  growth.  These  studies  of  the  effect  of 
layup  illustrate  a  major  difference  between  unidirectional  and 
multidirectional  laminate  delamination  behavior.  The  opportunity  for  more 
mechanisms  of  fracture  leads  to  a  greater  complexity  of  behavior.  It  is 
therefore  essential  that  attention  be  paid  to  the  fracture  morphology  when 
interpreting  the  data. 

Tests  were  done  to  investigate  the  dependence  of  J  on  crack  speed.  The 
results  for  a  unidirectional  layup  and  a  fiber-dominated  layup  are  given  in 
fig.  4.10.  In  both  cases  there  is  a  slight  decreasing  trend  in  J  with 
increasing  crack  speed.  This  trend  agrees  with  the  observation  that  when 
sudden  crack  jumping  occurred,  a  low  value  of  toughness  was  usually 
measured.  It  is  significant  to  note  that  no  strong  trend  can  be  established 
for  the  range  of  crack  speeds  experienced  by  one  specimen.  Thus  rate  effects 
are  not  important  in  interpreting  data  for  individual  specimens. 

The  importance  of  considering  the  fracture  morphology  when  interpreting 
the  delamination  toughness  data  has  already  been  mentioned.  Macro-  and 
microfractography  were  done  to  document  the  fracture  mechanisms  which  were 
observed.  Two  items  are  of  particular  interest.  First,  when  off-axis  plies 
were  present  at  the  fracture  plane,  toughness  values  were  sometimes  elevated 
even  in  the  absence  of  complex  fracture  mechanisms.  Apparently  this  was  due 
to  increased  surface  roughness  caused  when  bundles  of  fibers  were  pulled  away 
from  the  fracture  surface.  This  feature  was  only  observed  for  layups  with 
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off-axis  plies  at  the  delamination  plane.  For  a  more  brittle  system  the 
phenomenon  was  found  to  give  differences  of  a  factor  of  two  in  toughness 
between  unidirectional  and  off-axis  specimens  [ 4 . 7 j .  Figure  4.11  is  a  cross- 
section  of  a  specimen  (normal  to  the  direction  of  crack  propagation)  at  the 
fracture  surface.  A  bundle  of  fibers  is  separated  from  the  surface  at  this 
plane,  but  is  attached  to  it  at  another  point. 

A  second  phenomenon  which  is  of  significance  is  that  in  almost  all 
specimens  tested  delamination  was  found  to  occur  within  a  ply  rather  than 
between  two  plies.  Figure  4.12  shows  a  typical  situation  on  a  polished  cross 
section  at  a  fracture  surface.  Note  that  the  fracture  plane  does  not  pass 
through  the  interply  region,  but  stays  a  few  fiber  diameters  away  from  it. 
Practically,  the  presence  of  resin-rich  regions  and  ply  interfaces  has  only  an 
indirect  impact  on  the  mode  I  delamination  toughness  when  such  steady-state 
morphology  is  present. 

Additional  Studies:  Future  work  under  the  follow-on  grant  will  expand  on 
the  investigations  of  the  effect  of  geometry,  layup,  rate  effects,  and 
fracture  morphology.  Further  experiments  will  be  performed  to  explore  the 
effect  of  specimen  width,  especially  as  it  relates  to  crack  front  curvature. 
The  effect  of  the  number  of  plies  and  the  bending  stiffness  of  the  legs  of  the 
DCB  specimen  need  to  be  addressed.  Also,  tests  will  be  conducted  with  various 
fiber  angles  at  the  delamination  plane.  The  dependence  of  delamination 
toughness  on  crack  speed  for  matrix-dominated  layups  will  be  investigated.  An 
attempt  will  be  made  to  systematize  the  correlation  between  fracture 
morphology  and  delamination  toughness.  In  addition  to  these  continuing 
activities,  new  study  will  be  begun  to  understand  the  effect  of  continuum 
damage  on  the  apparent  toughness.  This  will  probably  include  delamination 
tests  which  simulate  mechanical  states  in  real  structures,  such  as  tests  using 


Figure  4.11 


Micrograph  of  a  polished  section  of  a  failed  DCB  specimen 
showing  the  fracture  surface.  A  bundle  of  fibers  has  been 
pulled  away  from  the  surface.  (150X) 


Figure  4.12  Micrograph  of  a  polished  section  of  a  failed  DCB  specimen 

showing  the  fracture  surface.  The  large  circles  at  the  top  are 
the  material  used  to  mount  the  specimen  for  polishing.  (100X) 


a  modified  DC8  specimen  with  bending  and  stretching  inputs.  Some  testing  of 
specimens  with  fatigue  predamage  is  planned.  In  these  studies  of  damage, 
emphasis  will  be  on  exploring  the  general  applicability  of  the  J  integral. 
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An  approach  using  strain  energy-like  potentials  to  char¬ 
acterize  deformation  and  fracture  of  Inelastic,  nonlinear 
composite  materials  Is  described  The  Inelasticity  may  be 
due  to  various  causes.  Including  microcracking,  microslip¬ 
ping.  and  rate  processes  responsible  for  fading  memory 
(viscoelasticity)  The  concept  of  work  potentials  is  intro¬ 
duced  first,  and  then  arguments  are  given  for  their  exist¬ 
ence  for  inelastic  materials.  Emphasis  in  the  paper  Is  on 
elastic  composite  materials  with  changing  or  constant 
states  of  distributed  damage.  Experimental  results  on  poly¬ 
meric  composites  are  subsequently  presented  to  illustrate 
this  approach  to  deformation  and  fracture  characteriza¬ 
tion  Finally,  extension  to  viscoelastic  behavior  Is  dis¬ 
cussed. 
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INTRODUCTION 

any  important  results  on  the  deformation 
and  fracture  of  linear  and  nonlinear  elastic 
materials  have  been  obtained  by  using  strain 
energy  functions  or  potentials  to  characterize 
material  response.  The  thermodynamics  of  re¬ 
versible  processes  provides  theoretical  support 
for  the  existence  of  these  potentials  and  iden¬ 
tifies  them  as  free  energy  and  interna!  energy 
for  Isothermal  and  adiabatic  processes,  respec¬ 
tively  (e.g..  Fung  (1)).  Besides  serving  as  the 
basis  for  powerful  methods  of  exact  and  ap¬ 
proximate  structural  analysis,  strain  energy 
functions  have  been  used  in  the  prediction  of 
effective  or  average  constitutive  properties  (or 
their  upper  and  lower  bounds)  of  linear  multi¬ 
phase  media  in  terms  of  properties  and  geome¬ 
try  of  the  phases,  as  reviewed  by  Hashin  (2). 
Included  in  the  many  publications  In  this  area 
are  studies  of  the  influence  of  small  distributed 
cracks  on  the  effective  stress-strain  behavior  of 
monolithic  and  composite  materials,  like  those 
described  by  Hashin  (2)  and  Kachanov  (3) 
Methods  of  characterization  and  analysis  us¬ 
ing  local  and  global  strain  energy-like  potentials 
for  certain  Inelastic  materials,  namely  viscous, 
plastic,  and  elementary  types  of  viscoelastic 
bodies,  have  been  discussed  in  an  early  work 
by  Hill  (4).  Constitutive  equations  normally  em¬ 
ployed  for  linear  and  nonlinear  viscous  bodies 
are  fully  analogous  to  those  for  elastic  media, 
in  which  strain  rate  replaces  strain  (4).  For  the 


linear  viscous  and  viscoelastic  cases  one  may 
use  irreversible  thermodynamics  (5)  or  special 
types  of  material  symmetry,  l.e..  cubic  and  iso¬ 
tropic  (6).  to  argue  for  the  existence  of  strain 
energv-like  constitutive  potentials  in  terms  of 
physical  or  Laplace-transformed  variables  (7). 
While  experimental  data  on  multiaxial  nonlin¬ 
ear  viscous  behavior  of  metals  (corresponding 
to  the  secondary  creep  stage)  may  be  character¬ 
ized  analytically  through  a  potential,  a  general 
theoretical  basis  for  this  constitutive  potential 
does  not  appear  to  exist  in  either  irreversible 
thermodynamics  or  in  models  of  micromechan¬ 
isms.  Rice  (8)  concluded  that  there  is  no  suffi¬ 
ciently  general  physical  model  of  slip  that  is 
capable  of  providing  a  firm  basis  for  the  exist¬ 
ence  of  a  creep  potential  Duva  and  Hutchinson 
(9)  give  a  good  illustration  on  the  use  of  poten¬ 
tials  to  construct  approximate  effective  consti¬ 
tutive  equations  of  nonlinear  viscous  compos¬ 
ites  In  this  analysis  the  composite  is  a  homo¬ 
geneous,  isotropic,  incompressible,  power-law 
nonlinear  material  with  a  given  dilute  concen¬ 
tration  and  size  of  spherical  voids  or  penny¬ 
shaped  cracks. 

In  using  a  potential  to  characterize  constitu¬ 
tive  equations,  it  is  often  sufficient  to  account 
explicitly  for  only  a  dependence  of  the  potential 
on  the  stress  or  strain  (or  strain  rate)  tensor.  If 
the  effect  of  different  temperatures  or  other 
parameters,  such  as  microvoid  or  microcrack 
fractions  and  sizes,  are  of  interest,  then  one 
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would  of  course  have  to  consider  how  these 
quantities  affect  the  constitutive  potential.  An 
example  of  such  a  potential  for  an  elastic  ma¬ 
terial  with  damage  is  the  volume-averaged 
strain  energy  density  of  a  material  specimen, 
D,).  where  ev  are  components  of  a  suita¬ 
bly  defined  volume-averaged  strain  tensor,  and 
Dk  represents  a  set  of  "damage"  parameters  that 
define  the  current  damage  state  (e.g..  micro- 
crack  sizes).  The  stresses  for  this  material  are 
then  obtained  by  differentiating  W  with  Dk 
fixed. 

su  ■  dW/dev  (1) 

In  Refs.  2.  3.  and  9,  the  effects  due  to  specified 
sizes  and  concentrations.  Dk.  of  microcracks 
are  considered.  If  It  Is  further  desired  to  char¬ 
acterize  the  effective  constitutive  behavior 
when  the  D»  change  with  time  as  a  result  of 
straining,  then  relationships  governing  these 
changes  must  be  determined.  For  discussion 
purposes  suppose  that  these  relationships  are 
known  and  that  one  may  solve  the  equations  so 
as  to  express  the  damage  parameters  Dk  In 
terms  of  the  instantaneous  strains  eu.  For  some 
cases.  It  may  then  be  possible  to  find  a  strain 
energy-like  potential  W[eu).  say.  from  which  the 
effective  constitutive  equations  can  be  derived 
by  differentiation, 

su  =  dW/de  y  (2) 

This  consitutlve  potential  would  depend  on  only 
the  instantaneous  strains  but  yet  account  for 
changing  damage  If  such  a  potential  could  be 
found.  It  would  be  like  that  used  to  characterize 
elasto-plastic  behavior  of  metals  by  the  Hencky 
deformation  theory  (10).  Similarly,  it  would  be 
analogous  to  the  potential  for  metals  discussed 
by  Rice  (8)  for  stationary  creep.  In  his  case  the 
"damage"  is  an  idealized  set  of  Internal  slips 
that  contribute  to  the  average  strain  rate  but  do 
not  appear  explicitly  in  the  effective  stress- 
strain  rate  equations  for  the  metal. 

The  present  paper  deals  In  large  part  with  the 
question  of  whether  potentials  analogous  to 
W[e{J)  exist  for  elastic,  viscous,  and  viscoelastic 
composites  with  changing  damage  (or,  more 
generally,  changing  microstructure);  emphasis 
is  on  elastic  behavior  with  damage.  Theory  and 
related  experimental  work  using  data  on  a  par¬ 
ticle-filled  rubber  and  fiber-reinforced  plastics 
are  discussed. 

We  should  add  that  there  are  already  many 
publications  on  the  thermodynamic  and  micro- 
mechanistic  bases  for  constitutive  potentials 
for  different  types  of  Inelastic  materials;  e.g.. 
Rice  (8.  11).  Coleman  and  Gurtln  (12).  and 
Schapery  (13  to  15)  However,  these  potentials 
depend  explicitly  on  "internal"  parameters 
which  reflect  the  microstructure  state,  and  are 
thus  like  W(e {l.  Dk)  In  Eq  1.  They  do  not  neces¬ 
sarily  lead  to  the  simplified  form  W(ey)  of  par¬ 
ticular  Interest  here. 


Strain  energv  potentials  have  been  used 
widely  in  fracture  mechanics  (e  fi  .  Broek  (16); 
For  fully  elastic  materials,  the  mechanical  work 
available  at  a  crack  tip  for  producing  an  incre¬ 
ment  of  crack  growth  is  equal  to  the  decrease 
In  potential  energy  (consisting  of  global  strain 
energy  and  the  boundary-work  potential)  Use 
of  this  relationship  has  resulted  In  remarkably 
successful  investigations  of  fracture  of  rubber 
In  Its  nonlinear  range  of  behavior,  which  are 
reviewed  by  Lake  (17).  as  well  as  fracture  of 
linear  elastic  materials  (16).  Andrews  (18)  as¬ 
sumed  a  strain  energy-llke  potential  exists  for 
rubber  with  hysteresis,  and  suggested  how  the 
hysteresis  would  affect  crack  growth.  When  a 
potential  exists  it  is  often  possible  to  use  Rice  s 
J  integral  theory  (19)  to  simplify  fracture  anal¬ 
ysis.  Schapery  (20)  recently  extended  the  poten¬ 
tial  energy  and  J  Integral  theories  to  elastic  and 
viscoelastic  materials  with  damage. 

Concepts  from  fracture  mechanics  are  used 
In  the  section  on  Multidimensional  Theory  to 
obtain  the  equations  needed  to  predict  micro- 
crack  growth,  and  thus  help  provide  the  basis 
for  potentials,  such  as  W(eu)  in  Eq  2.  Also, 
potentials  are  used  in  the  section  on  Delami¬ 
nation  in  Double  Cantilevered  Beam  to  account 
for  the  effect  of  inelastic  material  behavior 
(which  may  be  due  to  microcracking)  on  the 
growth  of  a  macrocrack  in  the  form  of  a  delam- 
inatlon. 

In  most  of  this  paper  it  is  assumed  the  mate¬ 
rials  are  elastic  when  the  damage  is  constant. 
In  the  Viscoelastic  Behavior  section,  a  special 
representation  of  viscoelastic  behavior  pro¬ 
posed  by  Schapery  (1 5.  20)  is  used  to  extend  the 
elastic  theory  with  damage  to  linear  and  nonlin¬ 
ear  viscoelasticity;  viscous  behavior  appears  as 
a  special  case. 

Finally,  it  should  be  mentioned  that  for  lack 
of  a  better  term  we  are  using  "damage"  when 
referring  to  characteristics  of  the  microstruc¬ 
ture  or  fabric  of  a  material  which  affect  consti¬ 
tutive  behavior  but  are  not  accounted  for  in 
elastic  or  fading-memory  viscoelastic  charac¬ 
terizations  of  continua.  Furthermore,  a  damag¬ 
ing  process,  such  as  considered  here,  could  be 
associated  with  crack  growth,  crack  healing, 
dislocation  creation  and  motion,  breaking  or 
reforming  entanglement  points  along  polymer 
chains  in  rubber,  etc.,  and  therefore  may  be 
structurally  detrimental  or  beneficial. 

ONE-DIMENSIONAL  THEORY 

The  definition  of  a  potential  for  elastic  mate¬ 
rials  with  damage  may  be  explicitly  introduced 
through  the  uniaxial  stress-strain  curve  in  Fig. 
1.  Let  us  suppose  that  a  previously  unloaded 
specimen  is  strained  monotonically  until  the 
strain  is  em.  (By  definition,  the  initial  state  is 
"undamaged.")  The  strain  Is  then  reduced,  as 
shown  in  Fig.  1 .  Assuming  that  the  bar  is  elast! 
and  has  constant  damage  during  the  unloading 
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The  strain  energv-like  quantity  W  is  acluailv 
,  the  net  work  to  the  materia!  at  any  suc<  of 

S  \  loading  or  unloading,  and  thus  it  will  be  called 

a  u'ork  potential  It  becomes  the  usual  strain 
/ /  energy  density  when  the  loading  and  unloadinc 

yj  j  /  /  curves  are  Identical.  Obviously,  a  work  poter.- 

,  -yf  I  //  /  tial  W  can  always  be  constructed,  given  the 

yf  /// //  uniaxial  stress-strain  behavior.  Eqs  3  and  4. 

\  yy  /  /  // //  Derivatives  of  multidimensional  equations 

\c  /  /  are  needed  in  the  next  section.  The  one-dimen- 

/y  /  /  /  sional  model.  Eqs  3  and  4.  is  useful  for  clarify- 

/  /  /  /  /  ing  some  of  the  analysis  ahead  of  time.  In  par- 

/  /  /  /  //  ticular,  observe  that  the  slope  of  the  unloading 

III  y\  curve  is. 

/  / If  )  em  dsv  dj 

/j  sU  de  de  de2  (  1 

/  / 1  m  The  loading  curve  s1  is  a  function  of  only  e. 

g  However,  when  using  the  upper  end  of  the  un¬ 

loading  curve  to  define  s^.  both  arguments  in 
Fig  1  Uniaxial  stress-strain  curve/or  material  with  In  j{£  ,  musl  ^  consldere(j  in  computing  the 

creasing  damage  during  loading  and  constant  damage  .  '  r  ° 

during  unloading.  derivative. 

+  (I0] 

stantaneous  stress  sL  and  using  de  de  dem  de  de  dem 

zation  as  Gurtin  and  Francis  (21)  ,  .  .  , 

naximum  strain  em  serves  to  de-  where  we  have  used  the  fact  that  em  -  e  on 

it  and  effect  of  damage.  the  loading  curve.  (The  last  term  in  Eq  10  is  the 

difference  in  slopes  of  the  two  curves  at  em  =  e). 
sL  =J{e.  em )  (3)  We  maV  rewrite  Eq  10  using  Eq  3  and  the 

second  expression  in  Eq  7  to  obtain 

J  curve,  the  maximum  strain  is 

•ain.  Hence,  viewing  the  loading  ds^  52U't 

it  at  the  upper  end  of  an  unload-  ~  ~de^~  +  dede~  I11) 

nay  write 

where  the  derivatives  of  are  to  be  evaluated 
s*-  =J(e.  e)  (4)  at  em  =  e  after  differentiation. 

t rn!n  t  MULTIDIMENSIONAL  THEORY 

to  an  arbitrary  strain  is 

p,  p.  For  characterization  of  multiaxial  stress- 

I  gL  de’  _  I  e') de'  (51  strain  behavior,  or  for  other  responses  that  de- 

Jo  J  pend  on  more  than  one  independent  input,  a 

,  ,  .  .  ,  ,  work  potential  does  not  necessarily  exist  How- 

ne  denotes  a  dummy  variable  of  ever,  it  can  be  expected  to  exist  for  some  real- 

te  net  work  input  to  the  sample  Jstic  situations,  which  will  be  discussed  here, 

nng  unloading  is  the  shaded  area  For  t^e  sa^e  Qf  generality,  let  us  use  as  inde¬ 

pendent  inputs  the  generalized  displacements 

X*  U  ~  1.  2.  . .  .  J).  The  responses  are  the  gen- 

svde'  eralized  forces  Q,,  which  are  defined  in  the 

(6)  usual  way  by  the  condition  that,  for  each  j. 

=  Wt(em)  +  J  J\e\  em)de'  -  QMj  (12) 

'■»  where  6(W7c)  is  the  virtual  work  input  associ¬ 
ating  loading  and  unloading,  re-  ?ted  wlth  vd?ual  displacement  bqr  Suppose. 

®  for  example,  that  we  let  each  qy  represent  a 

gradient.  dum/dxn  (m.  n  =  1.  2.  3)  of  a  three- 
V^/de.  su  =  dW^'/de  (7)  dimensional  displacement  field,  um.  and  let 

nt  to  let  W  denote  a  continu-  J™ k]  *  ^al  work  per  unit  initial  volume. 

vhich  equals  during  loading  Then  J  =  ,9'  afnd  f  <*  12  imP1,cs  the  Q-!*r\ 

equals  W*  during  unloading  components  of  a  stress  tensor  s  say  (for  large 

we  mav  write  for  both  loadinl  or  sma11  stralns'  FunS  I1))  In  order  to  charac- 

'  J  ®  terlze  the  behavior  of  laminates  using  classical 

P  plate  theory  (22).  one  may  want  to  associate  the 

s  *  dW/de  (8)  set  qj  with  the  middle  surface  curvatures  and 


period,  with  instantaneous  stress  sL'  and  using 
the  same  idealization  as  Gurtin  and  Francis(21) 
in  which  the  maximum  strain  em  serves  to  de¬ 
fine  the  amount  and  effect  of  damage. 

sL'  =/(e.  em)  (3) 

On  the  loading  curve,  the  maximum  strain  is 
the  current  strain.  Hence,  viewing  the  loading 
stress  as  a  point  at  the  upper  end  of  an  unload¬ 
ing  curve,  we  may  write 

s'-  =/(<?.  e)  (4) 

The  mechanical  work  (per  unit  initial  volume) 

during  loading  to  an  arbitrary  strain  is 

U*  =  U*(e)  =|  s1  de'  =  J  f(e\  e')de'  (5) 

where  the  prime  denotes  a  dummy  variable  of 
Integration.  The  net  work  input  to  the  sample 
at  any  time  during  unloading  is  the  shaded  area 
in  Fig.  1 . 


=  ^(e„ 


+  £  sude' 

J  J\e'.  em)de' 


Observe  that  during  loading  and  unloading,  re¬ 
spectively. 

s L  =  dV^/de.  sv  =  dW^’/de  (7) 

It  is  convenient  to  let  W  denote  a  continu¬ 
ous  quantity  which  equals  W1  during  loading 
(ew  =  e )  and  equals  W c  during  unloading 
(e  <  e„)  Then,  we  may  write  for  both  loading 
and  unloading  processes. 

s  -  dW/de  (8) 
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strains  In  this  case,  the  Q}  would  correspond  to 
moments  and  ln-plane  forces  per  unit  length, 
and  A(W'k)  would  be  virtual  work  per  unit  area. 

As  In  the  uniaxial  example,  we  assume  that 
when  the  damage  Is  constant,  the  body  (mate¬ 
rial  element,  test  specimen,  or  complete  struc¬ 
ture)  Is  elastic  in  the  usual  sense;  namely,  a 
work  potential  exists  with  the  property  that 

Qj-dVT/dqj.  (13) 

(Rather  than  using  the  terms  "loading"  and  "un¬ 
loading"  we  shall  Instead  now  refer  to  "damag¬ 
ing  processes"  and  "constant  damage  proc¬ 
esses."  since  we  do  not  w-ant  to  Imply  that  the 
damage  Is  always  constant  when  the  magnitude 
of  one  or  more  loads  or  displacements  decreases 
with  time.)  The  effect  of  damage  on  Qj  Is  as¬ 
sumed  to  be  fully  represented  by  a  set  of  "dam¬ 
age  parameters"  Fm  (n  ■  1 . 2. . . .  N)  in  the  next 
subsection. 

A  Special  Case 

Following  Schapery's  (20)  arguments.  It  will 
be  shown  that  for  a  suitably  chosen  W*  a  work 
potential  W°  exists  during  damaging  processes 
such  that 

Qj  =  dW°/d9j  (14) 


dQj  =  +  £  FdVV„  &F „ 

dqt  dqjdq,  <3F„  aqjdq, 

tf2  &wn  \  dFr,  arJ 

dFn2  +  dFndFcJ  dqj  dq , 

Clearly,  the  right-hand  side  of  Eq  18  is  the  same 
when  1  and  j  are  Interchanged,  and  therefore 
Eq  16  is  satisfied. 

Generalizations 

Extensions  of  VT  for  which  W°  exists  are 
discussed  in  Ref.  20.  and  W°  itself  is  given.  For 
example,  at  any  given  time  some  of  the  terms 
Wn  in  Eq  15  may  be  for  constant  damage  proc¬ 
esses  while  others  are  for  damaging  processes. 
Also,  the  potentials  may  depend  explicitly  on 
time,  and  thus  provide  for  effects  of  aging  or 
changing  physical  environments.  The  theory  in 
Ref.  20.  which  allows  for  large  deformations, 
uses  displacement  gradients  and  Plola  stresses 
instead  of  generalized  displacements  and 
forces.  However,  the  earlier  formulation,  in¬ 
cluding  its  extension  to  viscoelastic  behavior, 
carries  over  fully  in  terms  of  the  generalized 
variables  used  here.  Thus,  as  in  Eq  8.  we  may 
introduce  a  continuous  work  potential  W  for 
which 


where  W0  is  a  function  of  only  the  current 
values  of  One  special  W"  discussed  in  Ref. 
20  is. 

W1  =  Wj^)  +  I  W„(F„,  F„)  (15) 

n»  1 

where  W0  is  a  work  potential  without  damage 
effects  It  Is  assumed  that  all  of  the  functions 
F„  =  Fr(qj)  are  such  that  the  N  conditions  Fcn  = 
Fn  are  satisfied  simultaneously  during  all  dam¬ 
aging  processes.  (The  uniaxial  case.  Eqs  3  and 
4.  is  recovered  when  we  set  N  =  1  and  take 
U'„  =  0.  F,  *=  q ,  =  e.  Fc,  =  em.)  To  prove  that  W° 
exists.  It  is  necessary  and  sufficient  to  show 
that  the  generalized  forces  during  damaging 
processes  satisfy 


dQjfdq,  -  dQJdqj  (l J  -  1.  2.  ■  •  •  J)  (16) 

assuming  the  derivatives  in  Eq  16  are  continu¬ 
ous  (e  g  .  Greenberg  (23)).  The  forces  during  a 
damaging  process  are  taken  equal  to  those  in 
Eq  13  when  F„  «  F„  (which  Is  analogous  to 
saying  a  stress  on  the  loading  curve  in  Ftg.  1  is 
at  the  upper  end  of  an  unloading  curve).  Con¬ 
sequently.  we  may  evaluate  the  derivatives  in 
Eq  16  by  first  substituting  Eq  15  into  13, 


„  dW0  "  dWn  dF n 

q  m - +  y - 

dqj  dFn  dqj 


and  then  setting  F„  -  F„  and  differentiating  Eq 
1  7  (cf  Eq  11). 


Qj  =  dW/dqj  (19) 

even  if  the  damage  parameters  in  some  compo¬ 
nents  Wn  are  constant  while  others  vary  in 
time. 

Regardless  of  the  process. 
dQj/dq,  -  dQJdqj  -=  0 


except  at  the  points  of  change  from  one  process 
to  another,  considering  all  Wn.  The  derivatives 
in  Eq  20  are.  in  general,  discontinuous  at  these 
transition  points  (cf.  Ftg.  1  at  e  =  em)  and  thus 
Eq  20  does  not  apply  there.  Evidence  of  transi¬ 
tion  points  may  appear  in  experimental  data  as 
significant  but  somewhat  random  non-zero  val¬ 
ues  of  this  difference  of  derivatives  (for  !  #  j) 
over  short-time  Intervals.  This  experimental  be¬ 
havior  would  indicate  that  a  process  has 
changed  from  one  type  to  another,  damaging  to 
constant  damage  or  vice  versa. 

It  is  not  the  goal  here  to  develop  specific  phys¬ 
ical  models  that  gtve  rise  to  the  form  for  in 
Eq  15.  We  only  mention  that  one  based  on  a 
simple  microcracking  model  is  given  by  Schap¬ 
ery  (20).  Also,  for  characterizing  laminates  us¬ 
ing  classical  plate  theory,  each  Fcn  might  be 
proportional  to  a  ply  or  ply-pair  failure  surface 
(such  as  represented  by  the  Tsal-Wu  theory  (2). 
expressed  in  terms  of  the  local  ply  strains)  or 
other  local  invariant.  Since  the  local  strains  are 
(linear)  functions  of  the  mid-plane  strains  and 
curvatures,  one  obtains  F„  ■  F„(q,)  if  the  latter 
strains  and  curvatures  are  included  in  the  set 
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q,.  the  summation  In  Eq  15  would  extend  over 
all  plies 

It  is  not  necessary  for  the  constant -damage 
potential  to  have  the  form  In  Eq  15  for  the  work 
potential  in  Eq  19  to  exist.  For  example,  a  dif¬ 
ferent  form  for  W*  was  given  in  Ref.  20  which 
contains  the  Henky  deformation  theory  of  plas¬ 
ticity  (with  elastic  unloading)  Another  example 
is  given  in  the  next  subsection. 

Micro-  and  Macrocracking 

The  work  potential  W  in  £q  19  may  be  a 
constitutive  potential  in  the  sense  that  this 
equation  could  be  a  stress-strain  equation  for  a 
composite  or  monolithic  material.  Alternatively. 
W  may  be  the  total  work  input  to  a  structure 
under  a  general  set  of  boundary  displacements 
q,  whose  constitutive  response  is  defined  by  a 
work  potential  density.  In  either  case,  the  con¬ 
stitutive  potential  may  account  for  some  effects 
of  microcracking,  microvoiding.  slipping,  etc., 
through  the  damage  parameters  Fcn.  However, 
the  form  of  the  underlying  potential  for  con¬ 
stant  damage,  which  has  been  discussed  so  far. 
is  not  completely  general.  Also,  effects  of  mac¬ 
rocracking  (such  as  large-scale  delamination) 
have  not  yet  been  explicitly  Introduced.  Thus, 
it  is  of  interest  to  know  if  a  work  potential  exists 
when  there  is  macrocracking  and  a  relatively 
general  distribution  of  growing  microcracks. 
This  question  will  now  be  examined  by  embed¬ 
ding  additional  cracks  in  the  body  characterized 
by  U':  the  index  k  will  be  used  to  indentifv  each 
of  these  cracks,  assumed  to  be  K  in  number. 
The  cracks  may  have  a  wide  range  of  sizes,  but 
it  is  assumed  that  the  scale  of  the  crack  tip 
failure  process  zone,  which  determines  the 
work  required  for  increments  of  growth,  is  such 
that  the  local  material  surrounding  the  failure 
zone  can  be  approximated  as  a  continuum,  and 
that  the  effect  of  the  failure  zone  on  the  contin¬ 
uum  can  be  represented  by  tractions  acting 
along  the  local  crack  plane.  Then,  the  virtual 
work  equation  with  crack  growth  (see  Eq  13  in 
Ref  20),  which  applies  with  or  without  chang¬ 
ing  damage  in  the  continuum  and  regardless  of 
whether  or  not  growth  is  self-similar,  gives  the 
available  crack  tip  work  per  unit  of  new  surface 
area  as  -6W/6Ak-.  iW  is  the  change  in  the  work 
potential  for  the  total  body  due  to  the  Increase 
in  area  SAk  of  the  fc*h  crack  with  all  q,  fixed. 
Denoting  the  available  work  as  Ck,  we  may  thus 
wrtte 

Gk  =  -dW/dAk  (21) 

where  k  *  1 ,  2.  •  •  •  K;  also.  W  is  considered  to 
be  a  function  of  generalized  displacements  q, 
and  oriented  crack  areas  Ak.  The  quantity  Gk  is 
commonly  called  the  energy  release  rate.  The 
work  potential  may  also  depend  on  damage  pa¬ 
rameters  FCn  in  that  It  is  the  W  In  Eq  19  except 
for  the  fact  that  the  body  now  has  K  additional 
cracks;  the  virtual  work  (Eq  13  in  Ref.  20)  from 
which  Eqs  19  and  21  follow,  is  shown  in  Ref. 


20.  to  be  valid  with  crack  growth  in  bodies  with 
other  distributed  damaee 
In  order  to  predict  this  growth  we  aKf  need  t, > 
specify  the  work  required  Cr>.  say.  for  a  unit  of 
new  area  of  the  k'h  crack  area,  this  quannn  is 
the  so-called  critical  energy  release  rate  It  is 
not  necessary  to  assume  G(,  is  constant  or  is 
the  same  for  all  cracks.  However,  we  do  assume 
it  can  be  derived  from  a  fracture  work  potential 
VV}(  Ak).  where 

Gck  =  dV.'j/dAk  (22) 

For  example  if  the  critical  energy  release  rate 
for  all  cracks  is  constant,  but  not  necessarily 
the  same. 

K 

X  GckAk  (23) 

*-) 

Stable  (quasi-static)  growth  of  any  one  of  the 
K  cracks  is  specified  by  the  condition  that  re- 


quired  work  equals  available  work.  Gck 
thus. 

=  Gk; 

dWj/dAk  ~  —dW/dAk 

(24) 

If 

dWj/dAk  <  —dW/dAk 

(25) 

unstable  growth  occurs,  whereas  if 

aU>/(SAk  >  -dW’/dAk 

(26) 

there  is  not  growth. 

Returning  now  to  the  question  of  whether  a 
work  potential  exists  with  crack  growth,  we 
shall  see  that  it  does  for  the  model  defined  above 
if  the  growth  is  stable.  The  potential  is  denoted 
as  HV.  and  it  will  be  shown  that  it  is  simply  the 
work  of  fracture  plus  the  work  of  deformation 
of  the  elastic  or  Inelastic  continuum  W:  namely. 

Qj  «  dWT/dqj  (27) 

where 

WT=  Wj+  W  (28) 

The  proof  is  made  by  first  evaluating  the  deriv¬ 
atives  of  WV  while  allowing  for  the  stable  growth 
of  an  arbitrary  number  of  the  cracks;  hence 

wr  „  few,  aw\  dAk  aw 

dqj  ?  Ua*  +  dAk)  dq,  +  dqj  (  1 

For  those  cracks  that  do  not  grow.  dAk/Sq,  =  0. 
In  view  of  Eq  19  as  well  as  Eq  24  for  the  growing 
cracks,  Eq  29  reduces  to  Eq  27.  which  was  to 
be  shown. 

Equation  24  is  a  set  of  K'(0  <  K’  <  K)  equa¬ 
tions  whose  solution  gives  the  areas  Ak(q,)  of  K 
growing  cracks  in  terms  of  the  generalized  dis¬ 
placements.  If  we  assume  a  unique  solution 
exists  (at  least  for  small  changes  in  q, ).  then  the 
areas  (or  their  changes)  could  be  substituted 
into  Vtyand  W  in  Eq  28.  In  this  form  U>  would 
be  a  function  of  qj  and  only  those  Ak  that  are 
constant. 

Pursuing  this  representation  further,  let  us 

•7  - 
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suppose  for  simplicity  that  all  areas  Ak  simul¬ 
taneously  either  constant  or  vary;  denote  the 
constant  values  by  Ack.  Then  for  a  constant 
damage  process  (apart  from  possible  effects  of 
Fm  and  F„). 

W*T  =  WT(qj.  A**)  (30) 

and  for  a  damaging  process.  Ac*  — *  A*(qj). 

W?  =  Wr(q,.  Aklqj))  (31) 

Equation  30  does  not  necessarily  have  the  spe¬ 
cial  form  of  £q  15,  but  yet  a  work  potential 
exists  for  constant  and  changing  damage;  the 
limitation  is  instead  in  the  form  of  the  relation¬ 
ship  governing  A*(q,).  viz.  £q  24. 

Unstable  Crack  Growth 

Unstable  crack  growth  occurs  when  G*  >  Gc*. 
The  excess  work  predicted  from  quasi-static 
analysis  is  then  modified  by  dynamic  effects, 
and  the  quantity  WT  in  Eq  28  is  not  equal  to  the 
work  input  to  the  body.  This  does  not  necessar¬ 
ily  mean  a  work  potential  does  nto  exist.  In  fact, 
the  assumption  of  quasi-static  crack  growth 
was  not  used  to  arrive  at  Eq  19.  The  functions 
F„(q j)  may  reflect,  at  least  in  part,  an  average 
effect  of  unstable  rapid  steps  of  microcrack 
growth 

Significance  of  the  Areas  A* 

Through  principles  of  fracture  mechanics  we 
obtained  a  work  potential  WT,  and  furthermore 
related  it  to  physically  identifiable  parameters 
A*  and  material-related  functions  Wj  and  W. 
Conceptually,  all  cracks  are  considered  to  pre¬ 
exist  with  given  initial  sizes  and  orientations; 
but  many  or  all  may  be  so  small  initially  that 
they  have  no  real  effect  on  the  work  potential. 
In  principle,  as  many  A*  are  to  be  used  as  are 
needed  to  fully  define  the  instantaneous  loca¬ 
tion  and  orientation  of  all  crack  surfaces  and 
their  growth.  For  example,  a  crack  edge  that 
advances  nonuniformly  along  Its  length  may 
require  the  use  of  many  small  areas  or  param¬ 
eters  A*  to  define  the  changing  geometry.  With 
complex  arrays  of  cracks,  this  formulation  is 
not  practical  unless  idealizations,  such  as  peri¬ 
odicity  and  regular  shapes,  or  statistical  repre¬ 
sentations  are  Introduced.  Non-unique  solu¬ 
tions  A*(q j)  would  further  complicate  the  prob¬ 
lem.  giving  rise  to  effects  of  the  history  q,(f).  In 
such  a  case,  one  may  have  to  solve  for  changes 
In  A*  using  small  changes  in  q,.  Of  course,  our 
purpose  was  simply  to  argue  that  a  work  poten¬ 
tial  exists;  nevertheless.  It  should  be  recognized 
that  even  with  a  work  potential,  there  could  be 
effects  of  the  displacement  history.  Additional 
complications  could  arise  with  friction  between 
adjacent  crack  faces,  in  that  a  work  potential 
does  not  always  exist  if  there  is  appreciable 
energy  dissipation  through  sliding  processes; 
however.  It  should  be  recalled  that  significant 
plastic  deformation  (slip)  processes  may  occur 


in  metals  and  yet  a  potential  exists  for  some 
histories,  as  modeled  by  the  deformation  theory 
of  plasticity. 

One  could  think  of  the  parameters  A*  as  'in¬ 
ternal  variables."  such  as  used  in  irreversible 
thermodynamic  formulations;  they  need  not  be 
areas  as  long  as  they  fit  the  above  mathematical 
model.  Although  It  is  not  necessary  here,  we 
may  write  the  equations  that  govern  their 
growth,  Eq  24  or  equivalently  Gr*  *  G*.  in  a  rate 
form  similar  to  that  used  in  thermodynamic 
models.  First  differentiate  Gc*  *  G*  with  respect 
to  time. 


y  dGc*  dAm  y  dGk  dAm  y  dG *  dqj 

m  dAm  dt  '  t  dAm  dt  f  dqj  dt  ,J2) 

where,  for  simplicity,  explicit  time-dependence 
(e.g..  aging)  in  Gr*  and  G*  Is  omitted. 

Solve  for  dAm/dt. 


dAm 

dt 


1  Hmk 

kj 


dGk  dqj 
dq,  dt 


(33) 


where  H„*  is  the  Inverse  of  the  matrix  G*m 


G km  ■  dGCk/dAm  ~  dGk/dAm  (34) 


Complementary  Work  Potential 

In  many  problems  it  is  desirable  to  use  the  QJ 
instead  of  qj  as  the  independent  variables.  All 
of  the  preceding  theory  could  have  been  for¬ 
mulated  in  this  way.  in  which  a  complementary 
work  potential  Wc.  say.  would  be  used,  where 

qj  «  dWJdQj  (35) 

For  the  one-dimensional  case  In  Fig.  1.  Wc  is 
the  area  to  the  left  of  the  curves. 

Wc  -  J  eds'  (36) 


where  e  "'/(s',  sm)  ore* /(s',  s'),  depending  on 
the  curve  to  be  used,  and  sm  is  the  maximum 
stress.  The  relationship  between  WT  and  Wc  is 

(37) 
J- 1 

Observe  that  we  may  start  with  Eq  27  and 
then  define  Wc  by  Eq  37;  differentiation  of  the 
latter  equation  yields  Eq  35.  Alternatively,  we 
could  reverse  the  process.  Thus,  if  WT  exists  so 
does  Wc,  and  vice-versa.  It  should  be  noted  that 
these  potentials  are  multivalued,  and  therefore 
one  has  to  interpret  their  interrelationships  on 
a  process- by-process  basis.  For  instance,  in 
identifying  a  particular  (s.  e )  pair  for  the  ex¬ 
ample  In  Fig.  1,  It  is  obviously  necessary  to 
specify  whether  the  loading  or  unloading  curve 
is  to  be  used. 


ANGLE-PLY  COMPOSITE  BARS  UNDER 
AXIAL  AND  TORSIONAL  LOADING 

The  multidimensional  theory  provides  sup¬ 
port  for  using  work  potentials  to  characterize 
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the  behavior  of  materials  and  structures  with 
damage  In  this  and  In  the  next  two  Sections  we 
describe  studies  of  polymeric  composite  mate¬ 
rials  that  give  preliminary  experimental  confir¬ 
mation  of  this  characterization  for  damaging 
processes. 

Work  using  laminates  of  graphite  fiber-rein- 
forced  epoxy  under  axial  and  torsional  loading 
is  described  in  this  Section.  The  tests  specimens 
are  flat  bars,  as  illustrated  in  Fig.  2.  which 
consist  of  several  plies  or  layers,  each  being  a 
relatively  brittle,  elastic  composite  with  contin¬ 
uous.  unidirectional  fibers  having  an  angle  8 
with  respect  to  the  bar's  axis.  8  =  ±30°  are  used 
in  the  specimens  discussed  here.  The  unidirec¬ 
tional  material  was  supplied  in  pre-preg  form 
by  Hercules.  Inc.,  and  is  designated  as  AS4/ 
3502.  It  should  be  emphasized  that  even  though 
there  are  strain  gradients  and  consequent  non- 
uniform  damage  (primarily  in  the  form  of  dis¬ 
tributed  microcracks  and.  at  high  loads  edge 
delaminations)  the  multidimensional  theory 
may  be  used.  (Basic  stress-strain  behavior  using 
thin-walled  tubes  will  be  studied  in  the  near 
future  after  acquiring  additional  laboratory 
equipment.) 

The  generalized  variables  of  the  Section  on 
Multidimensional  Theory  will  be  identified  with 
the  specific  mechanical  variables  for  the  bar  as 
follows:  axial  elongation,  u  =  q,:  rotation  angle 
between  ends,  ft  =  q2:  axial  load.  F  =  Q and 
torque.  T  =  Q2.  The  total  work  potential  WT.  Eq 
28.  is  the  work  input  to  the  entire  bar  through 
the  relatively  rigid  grips 

The  necessary  and  sufficient  conditions  for 
existence  of  a  potential.  Eq  16.  for  the  present 
problem  reduce  to  the  single  equation. 

dF/dil  =  dT/du  (38) 

Before  using  Eq  38  with  experimental  data,  it 
is  helpful  to  replace  the  variables  by  measures 
of  stress  and  strain.  This  normalization  process 
eliminates  first-order  effects  of  specimen-to- 
specimen  dimension  differences.  Specifically, 
we  shall  use  -nominal-  stresses  and  strains  de¬ 
fined  as 

s  ■  F /he  t  ■  3 T/bc2  (3gj 

e  ■  u/L  7  ■  cft/L 

where  b  =  width,  c  =  thickness,  and  L  =  length 
(between  grips).  For  the  special  case  of  long, 
thin  homogeneous  specimens  (L»b»c)s  and 
e  are  the  uniform  axial  stress  and  strain  re¬ 
spectively.  and  r  and  7  are  the  in-plane  shear 
stress  and  strain  respectively  at  the  surface. 
This  is  shown  by  Timoshenko  and  Goodier  (24) 


F Ig  2  Laminate  specimen  used  In  axial-torsional  tests 


for  linear  isotropic  materials  for  orthotropic 
materials  whose  planes  of  symmetry  are  paral¬ 
lel  to  the  specimen  surfaces.  It  can  be  shown 
that  the  same  formulas  apply  except  the  width 
is  modified  by  a  ratio  of  moduli  The  shear 
strain  magnitude  is  zero  at  the  mid-plane  and 
increases  linearly  to  the  specimen  faces  ( L  x  b ) 
in  Fig.  2.  Whether  the  stated  conditions  apply, 
the  variables  in  Eq  39  are  useful  for  normaliz¬ 
ing  data.  Equation  38  becomes 

ds/dy  =  d(r/3)/de  (40) 

This  equation  has  been  used  to  analyze  the 
data  in  Figs.  3  and  4  by  first  writing 

t/3  =  To/3  +  g  (41) 

where  r0  =  1-0(7)  is  the  shear  stress  for  e  =  0: 
also  g  =  g[e.  7)  in  which  g(o.  7)  =  0.  Next, 
integration  of  £q  40  with  respect  to  7  yields 

d  fy 

s  =  —  g(e.  7  )dy’  +  sD  (42) 
oe 

where  s0  =  s0  (e)  is  the  axial  stress  when  7  =  0. 
Thus,  the  quantity 


Fig.  3.  Representative  shear  stress-strain  curves  with 
proportional  axial  and  torsional  straining  Hercules  AS4  ' 
35 02  graphtte/epoxy  24  piles  I30lt&.  2.5  in.  long  x  0.25 
In.  wide  x  0.13  In.  thick. 


Fig.  4.  Representative  axial  stress-strain  curies  with 
proportional  axial  and  torsional  straining  Jrom  the  same 
tests  used  for  Fig  3 


POLYMER  ENGINEERING  END  SCIENCE  MID-JANUARY,  198T,  Vo I  27.  Mo.  1 


M 


1^  !■  1  m  V  m  1  m  1  ■  V  m  'j  w  »'J  *'J  *\  w\*v: r: rjr 


i 

"  » 


<2 


a 


i 


* 

w\ 


K  A  Schapery 


d  l  y 

AS  “  de  X  9(e‘  y  )d>  (43) 

is  the  change  in  axial  stress  due  to  the  twist 
induced  shear  strain.  This  integrated  form  Is  to 
be  preferred  over  the  original  Eq  40  because  of 
the  inaccuracy  resulting  from  differentiation  of 
experimental  data,  especially  considering  the 
small  amount  of  data  available. 

The  procedure  used  to  check  for  the  existence 
of  a  work  potential  was  to  cross-plot  the  data  in 
Ftg.  3  so  as  to  obtain  g  (which  is  proportional 
to  the  change  in  shear  stress  due  to  axial  strain) 
as  a  function  of  y,  for  fixed  values  of  e.  and 
then  predict  the  modification  to  axial  stress.  Eq 
43.  Considering  the  limited  amount  of  data 
available,  it  is  desirable  to  curve  fit  analytical 
expressions  to  the  data  to  aid  the  needed  Inter¬ 
polations  and  extrapolations.  It  was  found  for  a 
wide  range  of  strains  that 

I  g(e.  y')dy' ~  AyaiSe,B*c-’'  (44) 
Jo 

where  A.  B.  and  C  are  constants.  Using  this 
expression  in  Eq  43  yields  the  change  in  axial 
stress  due  to  twist.  Only  for  e/y  ■=  0.92  is  there 
a  significant  effect  of  twist  prior  to  fracture;  the 
prediction  is  drawn  in  Fig.  4.  The  agreement 
between  theory  and  experiment  is  relatively 
good 

In  the  series  of  tests  shown  in  Figs.  3  and  4 
there  is  only  one  specimen  for  each  deformation 
history,  and  thus  the  small  differences  between 
most  curves  in  Ftg.  4  could  be  as  large  as  spec- 
Imen-to-specimen  differences.  Nevertheless.  It 
Is  encouraging  that  all  of  the  predictions  from 
Eq  43  turned  out  to  be  of  the  same  order  as  the 
observed  differences  In  axial  stress.  Axial  stress 
for  y  *=  0  Is  not  shown,  but  it  was  essentially 
the  same  as  for  e/y  =  1.78  in  Fig.  4  until 
premature  failure  occurred;  the  latter  results 
were  used  in  the  theoretical  predictions.  Al¬ 
though  not  needed  to  check  for  the  existence  of 
a  potential,  it  is  of  interest  to  observe  that  when 
there  Is  little  or  no  coupling  effect  of  twist  and 
axial  deformation,  the  stress-strain  curves  obey 
power  laws  over  a  wide  range  of  strains;  this  is 
shown  in  Ftg.  5  where 

s  -  em  and  r  —  yn  (45) 

We  have  conducted  additional  exploratory 
tests  using  proportional  straining  of  laminates 
with  various  widths,  thicknesses,  and  fiber  an¬ 
gles.  all  of  the  angle-ply  design  (±0)  with  bal¬ 
anced.  symmetric  layups.  The  behavior  is  sim¬ 
ilar  to  that  already  discussed,  with  comparable 
verification  of  Eq  43.  Close  to  the  end  point  of 
the  curves,  where  large  scale  delamination  or 
failure  at  the  grips  occurs,  theoretical  and  ex¬ 
perimental  curves  tend  to  separate,  as  seen  in 
Ftg.  4.  This  difference  may  be  due  to  inaccuracy 
in  the  extrapolations  needed  for  Eq  43  (consid¬ 
ering  the  small  number  of  specimens  used),  a 


change  from  damaging  to  constant  damage 
processes  and  vice  versa  (cf  discussion  of  Eq 
20],  or  an  inability  to  use  a  potential  In  a  highly 
damaged  state.  Future  studies  using  propor¬ 
tional  and  nonproportional  straining  should 
help  to  explain  this  behavior. 

For  some  layups  with  sufficient  twist,  mode 
m  edge  delamination  occurs  prior  to  significant 
material  fracture  near  the  grips.  As  a  result, 
properly  designed  bar  specimens  with  and  with¬ 
out  initial  delaminations  may  be  useful  for  stud¬ 
ies  of  this  type  of  delamination.  When  a  work 
potential  exists,  one  often  can  use  the  J  integral 
theory  or  energy  release  rates  to  account  for  the 
effect  of  distributed  damage  on  the  delamlna- 
tlon  growth. 

A  HIGHLY-FILLED  ELASTOMER  UNDER 
AXIAL  LOADING  AND  PRESSURE 

Several  years  ago  Farris  (25)  described  large 
deformation  studies  of  crossllnked  rubber  con¬ 
taining  65  volume  percent  Df  relatively  hard 
particles.  Specimens  in  the  shape  of  slender 
rectangular  bars  were  subjected  to  confining 
pressure  and  uniaxial  loading  He  used  revers¬ 
ible  thermodynamics  as  a  basis  for  predicting 
the  effect  of  pressure  on  the  axial  stress-strain 
behavior.  These  predictions  were  in  quite  good 
agreement  with  the  measurements,  in  spite  of 
strong  effects  from  the  irreversible  processes 
of  microcracking  and  void  growth.  Here  we  re¬ 
examine  the  behavior  and  use  the  present  work 
potential  theory  as  a  basis  for  making  similar 
predictions. 

The  specimen  and  relatively  rigid  grips  are 
shown  schematically  in  Ftg.  6.  This  assembly 
was  placed  In  a  chamber,  where  It  was  first 
pressurized  to  a  constant  value  p  and  then 
stretched  axially  at  a  constant  crosshead  rate; 
the  axial  force  acting  on  the  grips  Is  F.  Repre¬ 
sentative  stress-strain  and  dilatation-strain 
data  are  in  Fig.  7. 

In  order  to  select  the  generalized  variables  In 
Eq  12,  we  use  for  5(  Wk)  a  virtual  work  per  unit 
initial  specimen  volume. 


ooo  one  o  o 

hi  m*i» 


Fig.  5.  Logarithmic  plots  of  stress-strain  data  from  Figs. 
3  and  4  showing  power  law  behavior. 
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Fig  6  Particle  Jilted  elastomer  under  pressure  and  axial 
loading 
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Fig.  7  Stress-strain  and  dllatatlon-straln  behavior  of  a 
highly  filled  elastomer  165  volume  percent)  at  a  series  of 
hydrostatic  pressures  cf.  Farris  125). 


5(  UTc) 
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where  T  and  6u  are  the  surface  traction  and 
virtual  surface  displacement  vectors,  respec¬ 
tively.  and  V0  is  the  Initial  specimen  volume: 
the  integration  Is  over  the  instantaneous  area 
A  of  the  specimen  and  rigid  grips.  On  all  sur¬ 
faces  except  on  the  grip  ends  where  F  is  applied. 

T  =  -pn  (47) 

where  n  is  the  outer  unit  normal  to  the  surface. 
It  is  helpful  to  write  the  normal  traction  on  the 
grip  ends  in  the  form 

T„  -  T,  -  p  (48) 

where  7 ,  is  def ined  by  this  equation .  Integrating 
TndA  over  the  grip  ends  with  area  AB,  regardless 
of  whether  or  not  T„  is  uniformly  distributed, 
gives  the  axial  force  as 

F~  Fi-pA,  (49) 

where 

F,  -  T idA  (50) 

is  the  axial  force  above  that  due  to  the  pressure. 
Equations  47.  48.  and  50,  along  with  the  as¬ 
sumption  of  purely  axial  movement  of  the  grips, 
reduce  Eq  46  to 

F,  iV 

6{WK\  «  — ^  4u(  —  p  —  (51) 

V0  Vo 


where  6u,  and  iY  are  the  virtual  axial  elnnc.i- 
tion  and  volume  change  of  the  specimen  re¬ 
spectively  Let  us  now  choose  for  generalized 
displacements  the  nominal  axial  strain  c  and 
dilatation,  v.  defined  in  the  usual  way. 

q,  =  e  ■  u i/Lg.  q2  =  L'»V/V0  (52) 

where  u,  and  V'  are  the  increases  In  specimen 
length  and  volume  from  the  initial  unstrained 
state  (in  which  the  length  is  L0  and  volume  is 
V0).  Comparing  Eqs  12  and  51  we  see  that  the 
generalized  forces  are 

0i  «  s.  Q2  =  ~p  (53) 

where  s  ■  F i/A0.  the  "nominal  stress." 

As  in  the  previous  Section,  we  shall  use  Eq 
16  (with  t  «  2,  J  *  1)  in  Integrated  form  to 
determine  if  a  potential  exists.  Namely,  substi¬ 
tute  the  variables  from  Eqs  52  and  53  and 
integrate  with  respect  to  u. 


s(e.  u) 


-&  (I.  "*) 


+  s(e.  Do) 


where  o0  is  a  constant  reference  value  of  dila¬ 
tation.  This  result  is  equivalent  to  that  used  by 
Farris  except  for  an  additional  term  a.  .sing 
from  surface  or  fracture  energy,  which  he  at¬ 
tributed  to  the  formation  of  vacuoles.  However, 
he  subsequently  neglected  this  term  and  then 
used  the  theory  and  crossplots  constructed  f rom 
data  in  Fig.  7  to  make  the  predictions  in  Fig  8. 
For  s(e.  v0)  the  curve  for  p  =  500  psi  was  used: 
although  the  dilatation  is  not  constant  in  this 
case,  it  is  very  nearly  so  (cf.  Fig.  7). 

It  should  be  noted  that  Eq  54  is  correct  as  it 
stands,  in  that  surface  or  fracture  energy- 
changes  are  taken  into  account  implicitly  when 
one  uses  data  for  damaging  processes  The  un¬ 
derlying  potential  is  Wr.  Eq  28.  which  consists 
of  the  work  of  fracture  W }  plus  the  deformation 
work  W.  Furthermore,  no  restriction  has  been 
imposed  on  the  magnitude  of  the  strain  or  Its 
uniformity:  the  strains  are  in  fa.  t  quite  large, 
and  are  nonuniform  at  least  close  to  the  grips 


Fig  8.  Comparison  of  calculated  I points .  Eg  54)  and  mea 
sured  fllnesl  stress  strata  behavior  from  data  in  Fig  7  cf 
Farris  125). 
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The  agreement  between  results  from  experi¬ 
ment  and  potential  theory  ts  seen  in  Fig .  8  to 
be  quite  good.  The  discrepancy  is  relatively 
small  compared  to  the  differences  between  the 
various  stress-strain  curves  and  that  at  500  psi: 
the  error  that  does  exist  may  be  largely  due  to 
the  moderate  amount  of  viscoelasticity  exhib¬ 
ited  by  this  material.  Note  that  the  extent  of 
damage  is  large  at  the  low  pressures,  in  the 
sense  that  without  microcracking  and  the  sub¬ 
sequent  development  of  microcavities,  the  dila¬ 
tation  would  be  negligible  compared  to  the  val¬ 
ues  in  Fig.  7.  Also,  the  uniaxial  stress-strain 
curves  would  be  pressure-insensitive  since, 
with  Increasing  pressures,  the  curves  approach 
that  for  essentially  zero  dilatation. 

The  results  in  Figs.  7  and  8  are  from  tests 
conducted  at  constant  pressure.  Farris  also 
gave  the  results  in  Fig.  9.  which  include  a  test 
In  which  the  pressure  was  Initially  at  500  psi; 
while  the  sample  was  being  strained,  the  pres¬ 
sure  was  suddenly  lowered  to  40  psi;  following 
additional  straining,  it  was  increased  to  500  psi. 
It  is  seen  that  after  a  short  period  of  time  follow¬ 
ing  each  pressure  change,  the  stress-strain 
curve  tends  to  approach  the  one  for  which  the 
pressure  was  constant  during  the  entire  strain¬ 
ing  period.  In  other  words,  there  is  not  a  strong 
effect  of  pressure  history  in  this  case. 

DELAMINATION  IN  DOUBLE 
CANTILEVERED  BEAMS 

The  symmetric  split  beam  test  depicted  in  Fig. 
10  is  now  commonly  used  to  determine  the  crit¬ 
ical  energy  release  Cc  for  the  opening  mode  of 
delamlnatlon  of  fiber-reinforced  plastics.  When 
there  is  a  significant  volume  fraction  of  fibers 
that  are  not  parallel  to  the  beam  axis,  the  two 
legs  may  be  highly  inelastic,  thus  invalidating 
the  standard  elastic  methods  used  to  obtain  Gc 
from  experimental  data  on  load,  deflection,  and 
crack  length.  As  an  illustration  of  the  use  of  the 
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Fig  9  Stress  strain  curve  for  htghlyfllled  elastomer 
subjected  to  changes  in  the  pressure  levels.  After  Farris 
(25) 


F 


Fig.  10.  Double  cantilevered  beam  (DCB)for  delamlnatlon 
fracture  studies. 

potential  theory  in  the  section  on  Multidimen¬ 
sional  Theory  for  fracture  analysis,  we  shall 
derive  an  equation  for  determining  Gc  in  inelas¬ 
tic  beams. 

It  turns  out  that  the  complementary  work 
potential  Wc.  Eq  35,  is  used  more  conveniently 
for  this  problem  than  the  work  potential  In  the 
formulation  we  shall  employ,  the  potential  WT 
in  Eq  37  will  Include  the  fracture  work  of  all 
microcracks  but  not  that  of  the  delamination. 
For  this  case  the  work  that  becomes  available 
at  the  delamination  crack  tip  for  a  unit  of  new 
area  of  surface  A  (projected  onto  the  delamina¬ 
tion  plane),  is  given  by 

G  *  dWJdA  (55) 

where  th'  derivative  Is  taken  with  generalized 
forces  held  constant.  This  formula  can  be  de¬ 
rived  by  first  observing  that  the  total  variation 
of  Eq  37  may  be  written  in  the  form 


The  right  side  vanishes  by  virtue  of  Eqs  27  and 
35,  and  therefore  dWc/dA  ■  -dWT/dA.  Equation 
21  then  yields  Eq  55  since  we  may  use  Eq  21 
for  the  delamination  crack  in  which  the  work 
potential  is  WT  instead  of  W,  allowing  for  the 
work  of  microcrack  extension  introduced  in  Eq 
28. 

For  the  DCB  specimen  in  Fig.  10.  let  F  ~  Qt 
and  2u  «  q,.  Use  a  and  b  to  define  the  instan¬ 
taneous  crack  length  and  specimen  width  re¬ 
ferred  to  the  unstrained  (flat  beam)  geometry 
and  then  take  A  «  ab.  The  complementary  work 
for  the  total  test  specimen  is 

Wc  -  Wc(F.  a)  =  2  I  wcdxi  +  Wc'  (57) 
Jo 

where  xt  is  a  coordinate  along  the  specimen 
centerline  (which  defines  the  location  of  mate¬ 
rial  points  in  the  unstrained  geometry)  and  x  * 
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0  ts  at  the  Initial  line  of  load  application  The 
quantity  u\  is  defined  as  the  complementary 
work  density  (based  on  unit  length)  for  each  leg 
from  classical  beam  theory:  namely,  the  theory- 
based  on  the  assumption  that  material  planes 
that  are  initially  normal  to  the  beam  axis  re¬ 
main  plane  and  normal  when  the  beam  is  de¬ 
formed  The  correction  to  beam  theory  is  de¬ 
noted  by  Wc’.  The  energy  release  rate  becomes 


G 


l  dWc 
b  Sa 


l  “*a)  +  b 


2  f  • "  ew,  . 

bX  ~^dx' 


(58) 


after  neglecting  dWc'/da.  For  long  beams  Wc'  is 
primarily  from  distortion  of  the  beam  immedi¬ 
ately  to  the  right  of  the  delamination  crack  tip. 
and  dWc'/da  can  be  shown  to  vanish  if  the  beam 
is  long  enough  that  the  crack  tip  is  essentially 
isolated  from  the  ends. 

Equation  58  is  not  restricted  to  small  strains 
and  rotations.  As  such,  the  horizontal  distance 
a  '  in  Fig.  10  may  be  significantly  less  than  the 
actual  length  of  the  crack  surface  a.  However, 
we  shall  simplify  the  analysis  by  assuming 
small  strains  and  rotations  and  further  assume 
both  legs  are  symmetric  and  balanced  lami¬ 
nates  in  the  undamaged  and  damaged  states. 
With  these  conditions  there  is  no  mid-plane 
stretching  or  twisting,  and  beam  theory  gives 
the  simple  result 

u.y  =  f  kcLM'  (59) 

Jo 


where  k  =  fc(M')  is  the  curvature  as  a  function 
of  the  moment.  A  moment-curvature  diagram 
for  inelastic  material  would  be  similar  to  the 
stress-strain  diagram  in  Fig  1 ;  the  complemen¬ 
tary  work.  Eq  59.  is  thus  to  be  evaluated  taking 
into  account  the  multivalued  relationship,  as 
discussed  in  the  One-Dimensional  Theory  Sec¬ 
tion  The  local  beam  moment  is 


M  =  Fxi  (60) 

and  therefore  the  integral  in  Eq  58  vanishes  as 
we  is  independent  of  a.  Thus,  we  obtain 

G  *  ~  uyfa)  =  \  f  *  kiM]dM  (61) 
b  b  J  o 

where  Ma  =  Fa  is  the  crack  tip  moment. 

This  result  has  been  used  by  Jordan  (26)  to 
analyze  delamination  of  two  different  graphite/ 
epoxy  materia]  systems  with  plies  having  var¬ 
ious  fiber  orientations;  one  system  had  a  rub¬ 
ber-toughened  resin  and  the  other  a  brittle  resin 
system.  Figure  1 1  gives  a  typical  load-displace¬ 
ment  diagram  for  the  latter  one.  The  short  in¬ 
clined  lines  represent  periods  of  no  crack 
growth  following  sudden  Jumping  at  loads  along 
the  dashed  line.  Measurement  of  crack  length 
and  corresponding  loads  Indicated  that  the 
crack  tip  moment  at  the  beginning  of  each  jump 
In  crack  length  was  approximately  constant  for 
ail  tests  in  the  series.  Four-point  bend  tests 


were  used  lodevelop  the  moment -curvature  dia¬ 
gram  of  beams  with  the  layup  of  one-hai1  of  th< 
DCB  specimen,  t  e  ,  for  one  leg  Knowing  the 
moment  for  crack  growth  from  the  DCB  tests 
and  the  moment-curvature  diagram  from  the 
four-point  bend  tests  made  it  possible  to  derive 
the  critical  energy  release  rate  from  Eq  61 .  the 
average  moment  at  which  crack  jumping  oc¬ 
curred  for  each  DCB  specimen  was  used  for  Ma. 
This  calculation,  which  is  illustrated  in  Fig  12. 
gave  values  for  Cc  that  were  practically  the 
same  for  all  layups  of  each  of  the  two  graphite/ 
epoxy  systems.  Indeed,  the  Cc  for  unidirectional 
(0  =  0°)  laminates  was  close  to  that  for  the  layup 
with  multiple  fiber  orientations.  In  contrast, 
standard  data  analysis  based  on  beam  deflec¬ 
tion  and  load  gave  Cc  values  that  differed  con¬ 
siderably  for  the  several  layups;  some  multiple- 
fiber  angle  layups  had  apparent  Gc  values  over 
twice  that  for  0  =  0. 

These  findings  not  only  help  to  support  the 
underlying  potential  theory,  but  also  reveal  sur- 


orwtuc  DisrucDOJr.  2„  (in.) 

Fig  11.  A  typical  load-displacement  record  for  opening 
mode  loading  of  a  double  cantilevered  beam  with  off-axis 
piles  at  the  center  Interface.  Hercules  ASA/3502  graph- 
Ue/epoxy  After  Jordan  (261. 


K 

Fig.  12.  Moment-curvature  dlagramfor  loading  showing 
energy  release  rate.  C. 
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prisingly  simple  behavior,  considering  espe¬ 
cially  how  much  microcracking  develops  in  the 
specimens  whose  bending  Is  not  fiber-domi¬ 
nated.  The  layup  Insensitivity  of  Cc  seems  to 
Indicate  that  the  local  normal  interface  stress, 
rather  than  local  layup-induced  shear  stress,  is 
the  primary  factor  in  delamination.  Of  course, 
the  findings  are  from  a  limited  set  of  tests  on 
only  two  composite  systems,  and  therefore  one 
should  be  cautious  about  extrapolating  the  find¬ 
ings  to  other  laminates. 

Finally,  it  Is  of  Interest  to  consider  the  rela¬ 
tionship  of  the  above  results  to  Rice’s  J  integral 

(19) .  as  extended  to  crack  growth  in  inelastic 
media  with  large  deformations  by  Schapery 

(20) .  The  quantity  J  is  defined  by  a  contour 
integral,  which  becomes  for  the  beam  problem 
in  Fig.  13, 

J  •  X  -  (Ti  S; +  Ti  £f)dl]  1621 

where  w0  is  the  work  potential  density.  T(  and 
T3  are  tractions  along  C.  and  u,  and  ua  are 
displacements;  the  indices  Indicate  components 
in  the  X!  and  x3  directions.  This  equation  is 
valid  for  large  strains  and  rotations  as  long  as 
we  interpret  xt  and  x3  to  be  coordinates  of  the 
undeformed  geometry.  The  integration  is  coun¬ 
terclockwise  along  the  curve  C  in  Fig.  1 3.  which 
includes  top  and  bottom  beam  surfaces  and 
vertical  segments.  The  right  vertical  segment  is 
taken  far  enough  from  the  crack  tip  that  the 
material  is  unstressed,  and  thus  gives  no  con¬ 
tribution.  For  the  top  and  bottom  segments.  dx3 
*  T]  =  T3  =  0.  and  thus  the  only  contribution 
comes  from  the  left  segment.  Assuming  small 
strains  and  that  the  lef t  segment  is  close  enough 
to  the  crack  tip  that  we  can  use  small  rotation 
beam  theory,  yields. 

2  r  F  du3 

J«-Jo  WW-2  (63) 

where  M  and  du3/dxi  are  the  moment  and  slope, 
respectively,  at  the  left  vertical  segment.  Inte¬ 
grate  k  *=  d2u2/dx1 2  to  obtain  the  slope,  and  Eq 
63  reduces  to 

2  fM* 

J  -Z  \  kiM)dM  (64) 

D  Jo 

where  Ma  is  the  crack  tip  moment.  This  is  the 


Fig  13  Contour  used  to  evaluate  J  Integral. 


same  result  derived  by  Rice  ( 1 9)  for  a  split  beam 
under  end  moments  Ma.  It  is  seen  that  the  result 
is  Independent  of  the  location  of  the  left  inte¬ 
gration  segment  and.  in  fact,  is  that  in  Eq  63 
when  the  segment  is  located  at  the  crack  tip 
(where  du2/dxl  =  0).  It  should  be  mentioned, 
however,  to  obtain  this  path  independence  (i.e., 
derive  Eq  64  from  63)  It  was  necessary  to  as¬ 
sume  if  Jc(M)  is  multivalued  (cf.  Fig.  J)  that  the 
unloading  curve  is  the  same  for  all  left  vertical 
segments  used.  This  latter  condition  will  be  met 
for  all  material  (to  the  left  of  the  current  tip) 
that  had  experienced  the  same  maximum  mo¬ 
ment  when  the  crack  tip  passed  by.  Inasmuch 
as  the  experimental  results  discussed  above  in¬ 
dicate  the  maximum  moment  is  constant  (and 
recognizing  that  the  moment  decays  with  dis¬ 
tance  from  the  tip),  this  condition  is  met  all  of 
the  way  to  the  location  of  the  initial  crack  tip. 
To  the  left  of  the  initial  tip,  the  maximum  mo¬ 
ment  is  less  than  M„  and  one  finds  J  depends 
on  the  location  of  the  left  segment.  This  path 
dependence  in  beam  theory  is  fully  consistent 
with  that  predicted  from  the  exact  J  Integral  for 
a  continuum  with  variable  damage  in  the  re¬ 
gions  of  unloading  (20).  Some  unloading  may 
occur  in  the  continuum  very  close  to  the  crack 
tip.  possibly  causing  path  dependence  and  thus 
affecting  Eq  64.  Weatherby  (27)  used  a  finite 
element  analysis  to  study  this  dependence  for  a 
strongly  nonlinear  Isotropic  beam  and  found 
that  the  effect  on  J  is  negligible;  his  analysis 
predicted  a  J  value  very  close  to  that  In  Eq  64. 

Fora  crack  propagating  at  a  constant  moment 
in  a  long  laminate  that  is  initially  homogeneous 
In  the  Xi  direction,  the  state  of  stress  and  strain 
In  the  neighborhood  of  the  tip  Is  constant  in 
time.  This  is  a  type  of  “self -similar"  growth  and 
therefore  (19.  20) 

G  -  J  (65) 

Equations  61  and  64  agree  with  this  general 
result.  Equally  important.  Eq  64  was  derived 
without  assuming  small  rotations  (except  for 
the  neighborhood  of  the  tip).  Therefore,  in  view 
of  Eq  65,  we  may  conclude  that  the  formula  for 
work  available  at  the  crack  tip.  Eq  61,  is  valid 
even  when  the  DCB  legs  undergo  large  rotations 
and  possibly  high  axial  tensile  strains  from  the 
axial  component  of  F.  When  this  geometric  non¬ 
linearity  exists  one  should  use  a'  (cf.  Fig.  10) 
Instead  of  a  to  determine  the  crack  tip  moment. 

VISCOELASTIC  BEHAVIOR 

For  some  types  of  linear  and  nonlinear  vis¬ 
coelastic  materials  the  Multidimensional  The¬ 
ory  may  be  used  by  simply  replacing  the  gen¬ 
eralized  displacements  with  “pseudo  general¬ 
ized  displacements.*  qj!  where 

fl/“  ■  EjT1  f  E(t  -  r.  t)  &  dr  (66) 

*■'0  Or 

Here  qj  Is  the  physical  generalized  displacement 
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in  terms  of  the  dummy  time  variable  of  integra¬ 
tion  r  and  other  relevant  quantities  such  as 
coordinates  a  The  quantity  E  ~  E[t  -  r.  f)  is  a 
relaxation  modulus,  which  imparts  hereditary 
characteristics  to  the  deformation  behavior. 
The  second  f  argument  In  £  allows  for  aging 
and  other  kinds  of  time-dependence,  such  as 
may  be  due  to  transient  temperature  The  coef¬ 
ficient  Er  is  a  "reference  modulus."  which  Is  an 
arbitrarily  selected  constant  that  Is  Introduced 
so  that  qR  and  qj  have  the  same  units.  The 
basis  for  using  Eq  66  to  extend  deformation  and 
crack  growth  theory  to  viscoelastic  behavior  has 
been  given  by  Schapery  (15.  20).  Here  we  men¬ 
tion  only  that  it  is  an  exact  approach  for  linear 
isotropic  viscoelastic  materials  with  a  Poisson's 
ratio  that  is  constant  in  the  undamaged  state 
and  for  nonlinear  viscous  materials.  The  latter 
case  is  easily  deduced  by  noting  that  if  the 
relaxation  modulus  is  proportional  to  a  Dirac 
delta  function,  i.e.,  £  =  £r{*  4(t  -  r),  where  t,  is 
a  time  constant,  then  Eq  66  reduces  to 

qj *  =  tftdqj/dt  (67) 

The  pseudo  displacement  is  proportional  to  ve¬ 
locity  with  this  modulus  choice.  Replacement  of 
q,  by  this  qf  in  Eqs  19,  27.  and  35  converts 
them  to  equations  for  viscous  media. 

Elastic-viscoelastic  correspondence  princi¬ 
ples  were  developed  (15.  20).  which  then  lead  to 
the  development  of  the  theories  for  effective 
properties  of  composites  (15)  and  crack  and 
damage  growth  (20).  To  date  only  limited,  but 
encouraging,  experimental  work  on  nonlinear 
viscoelastic  materials  has  been  done  to  verify 
this  approach  based  on  pseudo  displacements 
(28.  29)  In  a  study  (29).  a  highly-filled  elastomer 
was  subjected  to  complex  uniaxial  loading  his¬ 
tories.  The  experimentally  measured  axial  dis¬ 
placement  u  of  tensile  coupons  was  converted 
to  nominal  strain  e  and  then  to  pseudo  strain 
eR  using  Eq  66.  and  the  axial  stress  was  plotted 
against  eR  for  constant  damage  states.  It  was 
found  that  this  method  of  plotting  data  essen¬ 
tially  eliminated  viscoelastic  effects,  thus  con¬ 
firming  the  approach.  For  damaging  processes, 
the  material  behaved  as  a  nonlinear  elastic  ma¬ 
terial  with  history-dependent  damage  when 
represented  in  terms  of  eR.  Guided  by  simple 
models  for  damage  and  some  test  data,  the  non¬ 
linear  viscoelastic  stress-strain  equation  was 
developed  and  then  verified  experimentally  us¬ 
ing  loading  histories  not  Included  in  the  char¬ 
acterization  process. 

Besides  employing  the  maximum  value  emR  to 
account  for  damage,  as  em  is  used  in  Eq  3.  a  so- 
called  Lebesque  norm  was  used  in  (29), 


where  I  |  denotes  absolute  value  and  q  is  a 
positive  constant;  the  damage  parameter  Lq 
arises  from  viscoelastic  crack  growth  theory.  If 


q  is  sufficiently  large  Lu  land  its  gener.il :,\r. .  i. 
to  multidimensional  problems;  ma\  ofum  It. id 
to  constitutive  equations  that  are  anal  v  .u*-  n 
those  discussed  in  Sections  on  One-Dimen¬ 
sional  Theory  and  Multidimensional  Theorv 
but  there  may  be  explicit  dependence  of  the 
potentials  W,  Wj.  and  W(  on  time,  and  thus  the\ 
would  be  like  those  for  aging  elastic  materials 
with  damage.  For  example,  eh  >  0  and  q  =  6  5 
in  Ref.  29.  and  in  this  case. 

Lq  -  eRt>,q  when  deR/dt  >  0  (69j 

and 

Lq  —  emRtWq  when  deR/dt  =  0  (70) 

If  eR  reaches  a  maximum  at  t  =  tm  and  then 
decreases. 

Lq  =  emRtmUq  (71) 

Behavior  like  that  in  Eqs  69  and  70  reflects  the 
growth  of  microcracks  until  the  "driving  force" 
eR  falls  below  its  maximum  value. 

CONCLUSIONS 

An  approach  using  work  potentials  to  char¬ 
acterize  deformation  and  fracture  behavior  of 
inelastic  materials  has  been  described.  Some 
experimental  results  on  polymeric  composites 
were  presented  to  illustrate  it  and  give  a  prelim¬ 
inary  verification  of  the  theory.  If  the  use  of 
work  potentials  to  account  for  the  effect  of  dam¬ 
age  and  other  types  of  inelasticity  is  further 
substantiated  in  future  studies,  one  may  take 
advantage  of  the  simplifications  that  come  f rom 
this  approach  in  the  mathematical  modelling  of 
both  deformation  and  fracture  behavior  When 
pseudo  displacements.  Eq  66.  can  be  used  to 
extend  the  time-independent  characterization 
to  nonlinear  viscoelastic  behavior,  additional 
experimental  requirements  and  mathematical 
model  complexity  are  not  much  more  than  what 
is  needed  for  linear  viscoelastic  behavior. 
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NOMENCLATURE 

a  =  crack  length. 

A.  Ak  =  crack  areas. 

b  =  beam  width. 

e  =  nominal  strain  (axial  extension/ 

initial  length). 
e„  =  strain  tensor. 

£(t  -  t.  t)  =  relaxation  modulus 

F  -  force. 

Fn  =  damage  function 

Fcn  =  damage  parameter 

G,  Gk  *  energy  release  rates 
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Cc.  Cck  =  critical  energy  release  rales 
J  =  J  Integral 

k  ■  curvature. 

L  m  beam  length. 

Lq  “  Lebesque  norm. 

M  -  moment, 

p  -  pressure 

qj  *  generalized  displacement. 

Qi  -  generalized  force, 

s  ■  nominal  stress  (axial  stress/initial 

cross  sectional  area). 

Sy  *  stress  tensor. 

t  -  time. 

»  torque. 

-  displacement. 

-  dilatation. 

«  Increase  In  volume. 

■  Initial  volume. 

-  complementary  work  density. 

-  complementary  work  potential. 

WT  “  work  potentials. 

-  fracture  work  potential. 

*  Cartesian  coordinates. 

■  nominal  shear  strain. 

*  nominal  shear  stress. 

■  angle  of  twist. 
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diameter  of  1  m  or  more  The  sapwood  is  pale  reddish 
xellow  and  clearly  demarcated  from  the  heartwood. 
which  caries  from  dull  reddish  brown  lo  deep  red 
Gri'wih  rings  are  distinct  and  relatively  wide,  but  the 
summerwood  bands  are  narrow  The  grain  is  straight 
and  the  texture  coarse  and  even  The  wood  is  very 
light  and  soft  and  fairly  easy  to  work  with  hand  and 
machine  tools  It  dries  well  without  degradation  but 
is  easv  to  cleave.  The  heartwood  weathers  well  The 
wood  has  a  pleasant  scent  because  it  contains  about 
2re  of  hvdroxysugiresinol.  and  is  permeable  to  pre¬ 
servatives  The  timber  is  used  in  building  construc¬ 
tion.  panelling,  ceiling  boards,  joinery,  furniture, 
poles  piles,  vats  and  boxes  Species  with  similar 
properties  and  uses  are  taiwania  (Taiwania  cryp- 
tomertoidcs) .  Chinese  fir  or  sha  mu  ( Cuninghamta 
smensts  or  C  lanceolata)  and  Formosa  fir  (C. 
knntshtn  These  species  are  important  plantation 
trees  in  China  and  Taiwan 

si’Rict  ( Picca  j ezoensis ).  also  known  as 
yeddo  spruce,  yezomatsu  (Japan),  gamunbi  namu  or 
ka-munpi  (Korea)  and  yulin-sdng  (China),  grows  in 
northeast  Asia,  northern  Japan.  Korea  and  Sakhalin 
It  is  one  of  the  most  important  plantation  trees  in 
Japan  and  is  a  conifer  that  reaches  a  height  of  3(1  m 
and  a  diameter  of  1  m  The  sapwood  and  heartwood 
are  not  demarcated,  both  being  pale  beige  or  pale 
vellow  Fine  pitch  streaks  are  found  frequently.  The 
grain  is  straight  and  the  texture  fine  and  even.  Planed 
surfaces  are  lustrous  The  wood  is  verv  light,  soft  and 
resilient  lr  machines,  seasons  and  finishes  well.  The 
wood  splits  easily  but  has  excellent  weathering 
properties  I'ses  are  in  house  construction,  interior 
trim,  siding  joinerv.  musical  instruments,  poles, 
ship-  and  boatbuilding,  vehicle  bodies  and  airplanes 
and  for  pulp  and  fiberboard  Species  with  similar 
properties  and  uses  are  Ffokkaido  spruce  or  akay- 
ezomaisut/'  ftiehni).  Siberian  spruce  {P  obovata), 
Japanese  fir  or  todomatsu  ( Abies  sclialtnensts).  Tai¬ 
wan  hr  t A  ka* akamtt) .  Taiwan  spruce  (P.  mor- 
risonu  nij  euan  sun  ( P  asperata )  and  others  in 
China  tP  tnanshurica .  P  pungsanensts ,  P.  lonat- 
ensts.  P  ituranuaca  and  P  retroflex) 

See  ah--<  Timbers  of  Southeast  Asia.  Timber  Resources  of 
Ihe  XX  orlj.  XX  t>od  An  Overview 


BlbllOgrap  ll\ 

Bn'wn  XX  H  19~8  Timbers  of  the  World.  Vols  3-5;  Timbers 
of  Southern  Asia.  Timbers  of  South  East  Asia.  Timbers 
of  the  Philippines  and  Japan  Timber  Research  and 
Development  Association.  High  Wycombe 
Hirai  S  (ed  l  I9K0  Encyclopedia  of  Woods,  Vols  1-17. 

Kanae  Shobo.  Tokyo  (in  Japanese] 

Rendle  B  J  (ed  )  1970  World  Timbers.  Vol  3:  Asia  and 
Ausoaha  and  Sen  Zealand  Benn.  London 
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Time-Dependent  Fracture:  Continuum 
Aspects  of  Crack  Growth 

The  growth  of  microcracks  and  macrocracks  is  affec¬ 
ted  by  time-dependent  material  behavior,  regardless 
of  whether  it  is  limned  to  the  crack  tip  neighborhood 
or  exists  throughout  the  bodv  This  article  is  con¬ 
cerned  with  specific  ways  in  which  rheological  pro¬ 
perties  of  continua  affect  both  initiation  and 
continuation  of  crack  grow  th  under  quasi-static  con¬ 
ditions  (i.e  .  when  inertial  effects  due  to  straining 
can  be  neglected)  The  highly  damaged  and  failing 
material  at  the  tip  of  growing  cracks  is  not  explicitly 
modelled;  rather,  emphasis  is  on  the  use  of  mechanics 
and  properties  of  the  surrounding  continuum  to  pre¬ 
dict  local  crack  tip  deformation,  stresses  and  the 
mechanical  work  available  for  extending  cracks 
Deformation  fields  based  on  bounded  rather  than 
infinite  strains  at  crack  tips  are  used  in  order  to  retain 
important  local  strain  rate  effects.  The  analysis  is 
simplified  by  using  small  strain  theory  for  the 
continuum,  while  allowing  for  large  strains  in  the 
failing  material  at  crack  tips 

1.  Constitutive  Equations 

Outside  of  the  highly  damaged  and  failing  material 
at  crack  tips,  the  deformation  behavior  is  assumed 
to  be  characterized  by  a  nonlinear  viscoelastic  consti¬ 
tutive  equation  in  the  form  of  a  single  hereditary 
integral  for  the  strain  tensor: 

(’  be', 

F,  =  £r  Du  -  r.  M - dr  ( 1 1 

J„  e  T 

The  quantity  e'„  is  a  second-order  tensor  which 
depends  on  material  properties  and.  in  general,  is  a 
function  of  the  suess  tensor  c r„.  material  point  x  . 
and  time  l: 

e'„  =  e',{crki.x„.i)  (2) 

in  which  ov  =  <t4(j,.  t).  with  all  indices  taking  the 
values  1.  2.  3  The  coefficient  £R  is  a  free  constant 
which  will  be  termed  the  reference  modulus,  it  is  a 
useful  parameter  in  discussing  special  material 
behavior  and  introducing  dimensionless  variables 
When  e'  is  used  in  Eqn.  ( I ),  the  time  argument  is 
specified  as  the  variable  of  integration;  that  is.  t 
should  be  replaced  by  r  where  explicitly  shown  in 
Eqn  (2)  and  in  the  argument  of  the  stress  Tosimplifv 
notation,  the  arguments  of  stress  and  strain  will  not 
be  written  out  unless  required  for  clarity  For  all 
cases  it  will  be  assumed  that  e ,  =  r'  =  cr,  =0  when 
t  <  0  and  Of  r  —  r.  t)  =  0  w  hen  t  <  r  To  allow  for  the 
possibility  of  a  discontinuous  change  in  e'.  with  time 
at  r  =  0.  the  lower  integration  limit  in  Eqn  (1)  and 
succeeding  hereditary  integrals  should  be  interpreted 
as  0  unless  indicated  otherwise 

The  quantity  f'  will  be  called  pseudo  strain  Its 
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v'. 


evplict  dependence  on  x„,  in  Eqn  (2 1  accounts  w hen 
necesvtrv  lor  material  nonhomogeneitv .  t  is  intro- 
t i ■ . v  ed  to  allow  for  aping  and  time-dependent  residual 
strains  (such  as  those  due  to  thermal  expansion 
(Vhapctv  lghli)  1  he  function /)(r  -  t.  r)  is  a  creep 
compliance.  it  provides  creep  under  constant  stress 
a'  well  as  other  hereditary  effects  under  time-varying 
stress  in  both  aging  and  nonagme  materials.  The 
significance  of  e'  and  D  will  be  show  n  bv  considering 
some  special  cases 

First,  however,  it  will  be  useful  to  rewri’e  Eqns 
(1)  and  (2  i  Hy  expressing  stress  in  terms  of  strain 
hisiors  Supposing  that  the  inverses  exist.  Eqn  (2) 
becomes 

cr,,  =  (rjel .x„.  t)  (2) 

where  f  is  related  to  the  physical  strain  through  the 

iio.  else  ol  Fqn  (  1 .1 


■  I  n  f  ■■ 

=  l- h  PU  -  r.  n - dr 

i<r 


(2) 


in  which  F.  is  a  relaxation  modulus,  its  relationship 
to  /)  is  eiven  h\ 


/ hr-  r.n  —  Fit' 
a  r 


.TUlr  =  Hit  -  r, 


(5) 


where  r  o  and  Hu  -  r  l  is  the  Heaviside  step 
furt.tion  tie.  Hu  -  r  )  =  (i  and  I  for  r  <  r,  and  for 
i  ■  r  respectively ). 

•\  linear  viscoelastic  material  without  residual 
s’  t  esse  -  which  is  isotropic,  homogeneous  and  has  a 
.  t’st  ■  r i  I’oisson  s  ratio  i  is  characterized  by  Eqn 
i  i  ill  w  c  use 

r'  =  £  *'  (t  I  •+  liar  -  i<7,.5  ]  (hi 

where  iS  t*  the  Kronecker  delta  and  the  standard 
sy’ui;  i turn  com ention  is  followed  in  w hich  repeated 
indices  imply  summation  over  their  range  For  a 
uni.iv.il  stress  state.  o:,*(J  and  all  other  o,  =  0. 
fqn  ill  for  f ; ,  becomes 


=  DU  -  r,  t) - dr 

or 


(  I 


I  If  i»,  =  o  Hit  -  fn).  where  rn  a  0  and  o,  is  constant, 
f  qn  i  ~  l  reduces  to  eu  -  DU  -  t.,.  Do .  Since  fn  n 
is  cusi  imunlv  termed  the  creep  compliance,  this 
name  shall  be  used  for  D  throughout  this  article 
Similarly,  if  r,  =  f..W(  r  - /0)  for  a  uniaxial  stress 
s'atc.  where  f.  is  constant,  one  finds  from  Eqn 
| -*  i  tha'  the  relaxation  modulus,  a,,  r , .  is  equal  to 
lu  t  .  n  W  hen  the  second  argument  in  D  and  £ 
in  F.qns  t  1 1  and  (4)  is  dropped,  so  that  Pit  -  r)  and 
I.  if  r)  appear  in  Eqns  (1)  and  (4|,  respectively , 
1  '’..mo  vis. ..elastic  behavior  for  a  nonaging  material 
is  Jes.ribed 

I  he  mechanisms  which  may  require  the  aging  form 


to  be  used  for  D  and  £  (e  g  .  D  =  Du  -  r.rn  are 
not  limited  to  chemical  processes  Fot  example  tins 
form  accounts  for  the  effect  of  transient  temperatures 
on  the  creep  compliance  and  relaxation  modulus, 
and  Includes  the  familiar  thcrmorheologivallv  simple 
behavior  of  polymers  as  a  special  case  It  should  he 
noted  that  the  expression  D(t-  r|  ts  sometimes  used 
instead  of  Du  -  r .  t)  in  characterizing  viscoelastic 
behavior  of  an  aging  material.  Although  both  lorms 
are  equally  general,  the  latter  is  used  here  as  it  o  a 
more  convenient  notation  in  equations  which  govern 
crack  growth 

Allowing  now  for  nonlinear,  anisotropic  and  non- 
homogeneous  media,  it  can  he  seen  that  for  the 
special  case  of  a  constant  relaxation  modulus.  £  = 
£„.  Eqn.  (4)  reduces  to  e'„  =  f„.  Thus.  Eqn  (2) 
becomes  the  constitutive  equation  for  an  elastic 
material  (in  that  the  current  stress  depends  on  the 
current  strain,  not  past  values  of  strainl  An  equiv¬ 
alent  result  is  found  by  using  D  =  £„'  in  Eqn  I  I  ) 
V  iscous  behavior  resuits  by  using  £  =  r.  £„<’)( r  -  r| 
in  Eqn.  (4)  (where  6(t  -  r)  is  the  Dirac  delta  func¬ 
tion  and  q  is  a  time  constant),  or  by  setting  D  = 
(/,  £r  )■’((  -  r)  in  Eqn  (1)  In  this  case  the  pseudo 
strain  ts  found  to  be  proportional  to  the  strain  rate, 
that  is.  f'  =  ,i/dt.  and  thus  the  current  stress 
becomes  a  function  of  the  current  strain  rate:  Eqn 
( 1 )  takes  this  form  after  integrating  it  by  parts,  then 
differentiating  and  inverting  the  result 

Hereditary  integrals  will  he  used  throughout  this 
article,  and  consequentlv  it  is  desirable  to  introduce 
abbrev  iated  notation  for  them  Specifically  ,  for  anv 
function  f  of  time. 


{ Ddf)mE*  Du 


of 

r.  M  —  dr 

f>  r 

if 


{E  d'l  =  E r1  £(r  -  r.  i)  —  Jr  (*i 
(•  r 


Thus.  Eqns  (1)  and  (4)  become 

f  „  =  { D  de'. |.  f’  =  i£  dr 


<gi 


2.  Correspondence  Principle 

The  close  relationship  between  mechanical  stales  ot 
nonlinear  elastic  and  viscoelastic  media  with  station¬ 
ary  or  growing  cracks  is  given  in  this  section  It  is 
stated  in  the  form  of  a  so-called  correspondence 
principle,  and  serves  as  the  basts  for  the  development 
of  crack  growth  theory.  First,  let  us  introduce  a 
reference  elastic  solution  crR  .  r  R  .  uR  corresponding 
to  the  case  in  which  D  '  =  £  =  £ „  Thtv  solution  is 
defined  to  satisfy  the  field  equations. 


=  0 


si  144 


1*  V* 
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=H  —  *  — 

OX  ox 
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f R  =  f’UTR .  x„. M 

(13) 

correspondence  principle  will  be 
It  the  instantaneous  ceomelrv  ( mclud- 

me  rr  k  k'  i  is  the  vime  for  both  elastic  and  v  iscoelastic 
;i!'  [  cl  surface  traction  T  =  a  n  (where  n  is 
the  normal  \  ect or )  be  a  specified  function  of  time  and 
position  i « inch  \  anishes  w  hen  r  <  0)  on  all  surf  aces  5 
I  hen  'hi  nonlinear  viscoelastic  solution  based  on 
E:qn  l  1  i  is 

a  =rrB.  t,  ={Dde*}.  u,  =  {DduH\  (13) 

w  here  me  v  ariables  w  ith  superscript  R  r  re  defined  to 
s  i t of \  Iiljps  <  I < ' V—  ( 12)  and  the  traction  boundary 
ior  !  ■■■  7  =  (rn  n  on  all  surfaces. 

V  •••  ■  I .  first  observe  that  the  viscoelastic  stress 
ri  !  i  I 'i  meets  the  condition  T  =  a  n  on  5 
i "  h  .  ■:  includes  the  instantaneous  crack  faces)  and 
s  : : i - : i •. '  equilibrium  Eqn  (1(1)  The  pseudo  strain 
f  “  .'bt<  Eqn  ill)  (compatibility  condition),  and 
. iearlv  the  viscoelastic  strain  in  Eqn  (13)  satisfies 
the  same  equation  without  the  superscript  R  (since 
/>  i-  independent  of  jr  and  the  hereditary  integral  is 
■I"  •  -peiaipr  •  Thus,  except  for  consideration  of 

tl’-;a  .  .  n  .uts  due  to  crack  growth,  the  proposed 
I  o  .  solution  meets  the  compatibility 
icq  iirvtw.s  \\  ith  or  without  crack  growth,  relative 
■$»f\  •,vr.t  nt  between  crack  faces.  Au  .  is  the  dif- 
f.  K'  ’  'he  displacements  in  Eqn  l  13)  evaluated 
on  ad  i.ern  faces 


iu  =  1  DrfAuR)  (14) 

■•‘b.:-.  \iC  is  the  displacement  difference  in  the 
tile:-.-.,  elastic  problem  Since  the  instantaneous 

c. .  all  cracks  in  the  reference  problem  is 

spe  '  ■  be  the  same  as  in  the  actual  viscoelastic 
\i.  iv  correctlv  predicted  to  vanish  until  the 
'  w hen  a  crack  tip  reaches  am  given  physical 
I  ■  '  follows  from  the  fact  that  AuR  =  0  at 

I  o  aiion  w  hen  r  f,  (assuming  prior  crack- 
iiil  ..  .  'v  ‘imng  of  crack  faces  has  not  occurred), 
•v"  1  m  'on  implies  the  hereditarv  integral  in  Eqn 

i  !  4  i  v  atiohes  when  l  <  i, 

I  he  correspondence  principle  may  be  generalized 
to  allow  for  specification  of  displacement  L,  on  some 
oi  all  surfaces  (Schupery  1984)  In  this  case,  the 
specified  surface  displacement  in  the  elastic  problem 
is  (  “*  =■  ■,/  Jt  )  and.  as  in  Eqn  (13).  elastic  and 
vis.  "•.To';,  stresses  throughout  the  continuum  are 
equal  wiih  stationary  and  growing  cracks. 


3  ( it  neralized  1  Integral  and  the  Crack  Tip  Mode l 

the  J  integral.  widelv  used  for  fracture  analysis  of 
tine-independent  materials,  will  he  generalized  in 
tho  section  for  subsequent  use  with  the  nonlinear 
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viscoelastic  materials  described  hv  Ei|u  i  3  ■  In  t in¬ 
elastic  problem  the  important  path  iiulepemieme  of 
the  J  integral  follows  in  part  from  the  thermo¬ 
dynamically  based  result  that  a  potential  H  exists 
with  the  proper! v  that  <r.  =("'11  rr  1  he  classical 
deformation  theorv  of  plasticilv  for  the  loatlmc  ot 
metals  may  be  expressed  in  terms  of  a  potential 
analogous  to  the  strain  energy  density  U  .  and  there¬ 
fore  the  J  integral  is  often  useful  even  if  large-scale 
plastic  deformation  exists,  as  long  as  there  is  no 
significant  unloading  from  the  plastic  state  W  ilhotn 
excluding  unloading,  equations  for  crack  growth  in 
nonlinear  viscoelastic  media  may  be  developed  and 
expressed  in  terms  of  a  J  integral  in  many  cases  il  o 
in  Eqn  (3|  can  be  written  as 

a,,  =  t)4>  rtf'  ( 1 3 1 

It  has  been  shown  that  the  potential  T  exist-,  for 
materials  under  general  strain  histories  il  titev  ate  at 
least  elastic  under  sudden  straining  (Sclt.ipcty  lu>> :  i 
For  linear  viscous  media  one  may  invoke  Ons.tger  s 
principle  to  establish  Eqn  (15)  Although  thermo¬ 
dynamic  arguments  alone  apparently  cannot  be  used 
to  justify  Eqn  (15)  for  nonlinear  viscous  behavior, 
the  standard  equations  used  for  secondary  creep  of 
metals  mav  be  written  in  this  form  (Leckie  and  II, iy- 
hurst  1 9~4 1  In  antilogy  with  elasticitv  theory,  d'  is 
called  the  pseudo  strain  energy  density 

Before  developing  a  fracture  theory  wlm.h  uses 
Eqn  (15).  certain  simplifications  concerning  the 
crack  tip  must  be  made  The  idealized  cr.uk  tip 
geometry  and  a  local  orthogonal  set  of  Cartesian  axes 
(x.)  are  shown  in  Fig  1  In  the  unstrained  state  the 
crack  surfaces  in  the  neighborhood  of  the  tip  are 
assumed  to  he  planar  and  to  coincide  with  the  los.tl 
x |.x ,  plane .  vv  here  x  is  perpendicular  to  the  plane  of 
the  page,  and  the  x-  axis  is  embedded  in  the  con¬ 
tinuum  at  any  convenient  horizontal  location  The 
crack  tip  or  edge  is  a  straight  or  curved  line  in  sp.uc 
whose  intersection  with  the  plane  ol  the  page  is 
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Figure  I 

Cross  section  of  crack  in  nonlinear  v  iso  tetanic  m  venal 

showing  comours  ( - (  used  m  line  iricetalv  ( i"h  the 

opening  mode  of  displacement  is  drawn  ,  1 1 1 1  - .  >  i :  c  t  -  the 
basic  formulation  allows  (or  shearing  delomution 
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indicated  by  the  point  P  It  is  assumed  the  tip  is 
essentially  straight  and  parallel  to  x-,  in  the  neigh¬ 
borhood  of  P  (By  definition,  the  neighborhood  of  P 
refers  here  to  the  material  contained  in  a  sphere 
w  hieli  is  centered  at  P  and  has  a  radius  of  about  ten 
times  the  failure  zone  length  a  ) 

The  region  designated  as  the  failure  zone  in  Fig  1 
is  where  intense  damage  and  material  separation 
occurs  Outside  this  zone  it  is  assumed  that  Eqn.  (15) 
applies  This  term  is  preferred  over  other  common 
names,  such  as  process  zone  and  damage  zone, 
because  the  material  outside  may  be  damaged  or 
otherwise  changed  and  the  effects  taken  into  account 
in  a  more  general  theory  (Schapery  1984)  Also,  for 
analysis  purposes  it  is  more  convenient  to  define  P 
as  the  crack  tip  location,  instead  of  the  left  end  of 
the  failure  zone.  The  latter  position  is  called  the 
apparent  crack  tip 

Consider  now  the  line  integral 

£  *  <j!  ^ <l>  dx i  -  T,  ~~  ds  j  (16) 

where  the  path  integration  is  in  the  x,.t;  plane  and  is 
any  closed  contour  on  and  inside  of  which  <J>  exists 
and  which  does  not  enclose  cracks,  one  such  path  is 
that  designated  as  CT  in  Fig  1  This  integral  v  anishes 
if  (a)  a  stale  of  plane  strain,  plane  stress,  or  antiplane 
strain  exists  and  (b)  r» st>  Ax ,  =  f*<t>  Tt,  =  0  on  and 
inside  C'i  According  to  the  second  condition,  the 
most  general  admissible  form  of  <t>  is  <t>  =  <J>(f'. 
<  - .  / ):  dependence  of  r '  on  x ,  is  permitted,  but  any 
material  nonhomogeneity  is  restricted  to  variations 
normal  to  the  local  crack  plane  Note  that  J  =  (I 
follows  directly  from  elasticity  theory  for  the  con¬ 
tinuum  within  the  contour  CT  (Rice  19b8)  since  £  is 
expressed  in  terms  of  mechanical  state  variables  for 
the  reference  elastic  problem  and  the  potential  <t>  is 
defined  by  Eqn  (15);  recall  that  er*  =  a  and  fr  = 
e'  according  to  the  correspondence  principle. 

A  generalized  J  integral,  denoted  by  JK.  is  now 
introduced: 

J,  =  (  (<t >dx,  -  T,—ds)  (17) 

J c,  dx,  ' 

where  C,  is  the  portion  of  CT  indicated  in  Fig  1; 
integration  starts  at  point  1  and  is  taken  counter¬ 
clockwise  to  point  2.  In  contrast  to  the  early  elasticity 
theory  (Rice  1968).  it  is  not  assumed  that  the  crack 
faces  are  traction  free,  and  thus  that  part  of  C,  which 
is  adjacent  to  the  surfaces  is  retained  (A  later  study 
by  Palmer  and  Rice  (1973)  included  body  forces  and 
crack  surface  tractions  in  analyzing  slip  surfaces  in 
clay  without  time  dependence  )  The  condition  £  = 
0  rray  be  w ritten  as 

yv=y,  (18) 


=  f 


*uR 

T  - ds  . 

(IX, 


and  integration  starts  at  point  1  and  proceeds 
counterclockwise  along  the  edge  of  the  failure  zone 
to  point  2  The  path  C,  for  J,  is  arbitrarv  except  that 
it  starts  and  ends  at  the  apparent  crack  tip  ( points  1 
and  2)  It  is  not  necessary  to  use  a  contour  having 
segments  parallel  to  the  crack  faces;  but  if  the 
tractions  T,  vanish  on  such  segments  to  the  left 
of  the  failure  zone,  there  is  no  contribution  to  7, 
(since  dx:  =  T,  =  0)  and  therefore  C,  could  start  and 
end  anywhere  along  the  lower  and  upper  crack 
faces. 

Equation  (18)  provides  a  basic  relationship  be¬ 
tween  the  mechanical  state  of  the  continuum  through 
J >  and  the  characteristics  of  the  failing  material  at  the 
crack  tip  through  J,  Given  that  real  stresses  and 
strains  (including  those  at  P)  are  not  infinite,  in  many 
cases  one  can  neglect  the  contribution  to  J,  from  the 
tip  segment  encircling  P  and  assume  the  tractions  T 
are  equal  (with  opposite  signs)  across  the  upper  and 
lower  parts  of  C2;  such  simplicity  obviouslv  exists  if 
the  failure  zone  is  thin  (in  the  a;  direction)  relative 
to  a.  Then.  Eqn  (19)  becomes 


[“  I  dAu*  8An? 

'<  =  |  (r„ 


. 

+  r,y  — —  I  dl  1 2" ) 

where  o,  is  the  normal  stress  and  rM  and  r,,  are 
the  shearing  stresses  in  the  r,  and  o  directions, 
respectively,  along  the  interface  between  the  failure 
zone  and  continuum.  A uH  is  the  relative  displace¬ 
ment  vector  between  initially  adiacent  material 
points  on  the  crack  plane  The  normal  stress  and 
displacement  and  the  local  coordinate  £  =  a  -  x  ,  are 
indicated  in  Fig  2. 
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Figure  2 

Normal  stress  and  displacement  along  crack  plane 
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4  Crack  Growth 

Essential  features  of  equations  for  predicting  the 
time  at  which  crack  growth  initiates  and  subsequent 
crack  propagation  speed  a  will  be  illustrated  by 
assuming  a  slender  failure  zone  (compare  with  Eqn. 
(2111)  and  the  opening  mode  of  crack  tip  displace¬ 
ment  In  this  case,  sliding  displacements  along  the 
crack  plane  are  zero  in  the  neighborhood  of  the  tip 

4  I  Initiation  of  Growth 

In  order  to  avoid  excessive  mathematical  complexity 
in  illustrating  the  prediction  of  r,.  the  failure  stress 
distribution  o,  will  be  assumed  independent  of  £, 
and  denoted  by  om  Equation  (20)  becomes  7,  = 
o„,Ai/?3.  where  Au?,  is  the  elastic  opening  dis¬ 
placement  at  the  fixed  location  £=a:  the  tip  P. 
rather  than  the  apparent  tip.  may  move  when  r  «  (,  in 
our  crack  tip  model  Usings  =  J,.  the  displacement  is 
An?,  =  7,  om.  and  the  correspondence  principle, 
Eqn  (13).  yields  the  time-dependent  displacement, 

Am:.  =  {Dd(Jjom))  (21) 

where  o„.  is  not  necessarily  constant  in  time  if  the 
criterion  for  initiation  is  based  on  a  critical  value  of 
opening  displacement.  Ai/C.  it  is  possible  todetermine 
the  time  at  which  An;,  in  Eqn.  (21 )  is  equal  to  A i/c. 

An  alternative,  commonly  used  criterion  is  based 
on  mechanical  w  ork  In  this  case,  work  input  from  the 
continuum  is  equated  to  work  required  for  failure: 


viscous  material.  D  =  (/,£„)  '(/  -  r)  where  t.  o  a 
time  constant,  as  noted  previously  After  integrating 
Eqn.  (24)  by  parts,  it  becomes 


w  here  f  is  the  mechanical  work  per  unit  surface  area 
required  to  break  a  material  column,  with  cross- 
sectional  area  dA .  at  the  left  end  of  the  failure  zone, 
designated  by  the  subscript  a  on  the  stress  and  dis¬ 
placement.  a  factor  of  2  is  used  to  account  for  the 
two  units  of  surface  formed  for  each  unit  of  cross- 
sectional  area  of  the  column  The  so-called  fracture 
initiation  energv  T,  may  be  affected  by  the  history  of 
Air-.,  depending  on  the  characteristics  of  the  failure 
zone  For  the  case  in  which  o{  =  om,  Eqn.  (22) 
becomes 

f  om[d{D  d{J  Jo m))/dt]  dt  =  2T,  (23) 

If  is  independent  of  time  and  the  explicit  rep¬ 
resentation  for  the  hereditary  integral  is  used  (com¬ 
pare  with  Eqn.  (8|)  then 


D(t,  -  r.  r, 


dJ  s 

\  —  dx  =  2T, 
dr 


n/ijf' 


Jsdt  =  2E, 


It  may  be  observed  that  the  failure  zone  length  a  has 
not  been  assumed  small  in  Eqns  (21-25)  Also,  a 
generalization  of  Eqn  (21 )  is  easily  obtained  if  o,  = 
0i(Au?.r);  one  finds  that  a  function  of  7,  and  r 
replaces  7v/om. 

4.2  Crack  Propagation  Speed 

The  failure  criterion  based  on  mechanical  work  may 

be  written  in  the  form 


'»  5Air- 

o,  :  d:  = 

J  r\  ("*  £ 


For  an  elastic  continuum,  where  £>  =  £r'.  this 
equation  reduces  to  7,  =  2 E, .  the  failure  zone  may 
be  viscoelastic,  and  therefore  2I~,  is  not  necessarily 
constant  If  the  continuum  is  a  linear  or  nonlinear 


where  the  fracture  energy  T.  just  as  T,.  is  the  work 
required  per  unit  area  to  rupture  a  column  of  material 
with  cross-sectional  area  dA  in  the  failure  zone  In 
contrast  to  the  initiation  problem,  this  column  is  not 
initially  al  the  apparent  (free  surface)  crack  tip. 
instead,  the  column  is  initially  within  or  ahead  of  the 
current  failure  zone,  and  its  rupture  is  coincident 
with  the  arrival  of  the  apparent  tip  The  left  side  of 
Eqn  (2b)  times  dA  is  the  total  work  input  to  the 
column,  starting  with  the  time  the  tip  P  (see  Fig  1 1 
first  arrives  For  convenience  in  analyz.ng  continuous 
growth,  the  variable  of  integration  has  been  changed 
from  t  (see  Eqn  (22))  to  £.  where  ;  =  a  -  r,.  and 
the  differentiation  and  integration  in  Eqn  (2b)  are 
made  with  x,  held  constant  The  time-dependent 
opening  displacement  at  any  point  in  the  failure  zone 
is  A u-  =  {Dr/Au?}.  which  is  to  be  substituted  into 
Eqn  (26)  The  result  for  predicting  d  is  rather 
complex,  even  if  A u:  is  used  in  a  critical  displacement 
criterion,  because  o,  may  depend  on  £  or  directly  on 
Au;  or  An?,  and  o,  affects  the  displacement 

The  problem  can  be  greatly  simplified  if  (a)  o,  and 
At/?  for  fixed  £.  a  and  a  are  essentially  constant 
during  the  time  interval  a/d  when  the  crack  propa¬ 
gates  a  distance  o  and  (b)  aging  changes  in  the  creep 
compliance  during  this  same  interval  are  negligible 
With  these  conditions  and  the  observation  that  the 
relation  between  Ai/;  and  Au?  is  the  same  as  for  a 
linear  viscoelastic  material,  an  approximate  evalu¬ 
ation  of  the  hereditary  integral  may  be  used  (Schap- 
ery  1975  pp  369-88)  to  find 

At/,  =  £„D(r.  r)Au?  (2*) 

where  r  =  k£/a  and  k  is  a  factor  which  is  a  very  weak 
function  of  the  slope  n  *  d  log  D  d  log  r;  this  factor 
is  practically  1/3  for  the  entire  range  of  slopes 
(Osrjsl)  encountered  in  practice  For  a  linear 
viscous  body  (n  =  1).  k  =  1/3  exactly  if  £  =  a  and 
Of  =  om  Equation  (27)  stems  from  the  smooth  cusp- 
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shaped  opening  displacement  predicted  for  linear 
materials  and  the  typically  weak  dependence  of  n  on 
log  i  Assuming  the  strains  are  continuous  at  the  tip 
P.  the  helms  lot  of  the  local  displacement  in  nonlinear 
media  is  expected  to  he  close  to  that  for  linear  pre¬ 
stressed  media  (Brockway  and  Schapcry  1978)  and 
theiclore  cusp-shaped  Thus,  although  more  careful 
studs  of  this  behavior  is  needed  for  highls  nonlinear 
materials,  the  cusp  shape  is  expected  to  exist  in  many 
cases,  permitting  the  use  of  Eqn.  (27)  and  a  constant 
salue  for  k  Substituting  Eqn.  (27)  into  Eqn.  (26) 
and  using  the  same  type  of  approximation  as  in  the 
linear  theory  ( Schapery  1975  pp  369-88 ).  which  does 
not  require  o,  to  be  spacewise  uniform,  there  results 

[‘  <?A  u? 

EpD(kad.t)  a, — —  d:  =  2r  (28) 
Jo  d  = 

In  siess  o|  Eqns  (18)  and  (201.  the  integral  in  Eqn 
IT's l  mas  he  replaced  by  7,  to  derise  a  very  simple 
result  (sir  predicting  a: 

EpD(ka  a .  t)J >  -  2E  (29) 

The  integration  and  differentiation  in  Eqn.  (28) 
are  for  fixed  \ while  those  in  Eqn  (20)  are  for  fixed 
r  (i  e  .  fixed  crack  tip  location),  however,  for  the 
steady -state  propagation  assumed  during  the  time 
a  it.  these  integrals  are  equal.  Also,  r  in  the  second 
ai  gument  of  O  accounts  for  aging  in  the  creep  com- 
pliance.  and  the  aging  is  assumed  small  during  the 
time  inters  al  n  d  In  spite  of  the  assumed  constancs 
of  a  and  aging  during  a  a.  they  mas  vary  appreciably 
during  the  total  time  of  crack  propagation  Likewise, 
the  si?e  n  is  next  necessarily  constant,  as  it  is  related 
to  other  variables  such  as  7,.  its  prediction  is  illus¬ 
trated  in  Sect  5  for  a  nonlinear  power  law  material 
Abo.  it  has  not  been  necessary  to  assume  a  is  small 
in  derising  Eqns.  (27)  and  (29);  it  is  sufficient  to 
assume  the  change  in  0|  and  An?  is  small  due  to  the 
chance  in  distance  between  the  crack  tip  and  nearby 
gc  'metric  features  (such  as  another  crack  tip)  during 
the  generic  period  a  a 

The  similarity  of  the  formulas  for  initiation.  Eqn 
(24)  and  continuous  growth.  Eqn  (29),  is  note¬ 
worthy  As  a  special  case,  suppose  7,  is  constant  for 
t  ■  (l  Then  Eqn.  (24)  reduces  to 

E„D(r,.r,)71  =  2T,  (30) 

This  result  and  Eqn.  (29)  yield  l,  =  ka'a  if  T  =  T, 
and  aging  is  negligible  Thus,  with  k  =  1/3.  the  time 
required  for  initiation  is  only  one-third  of  that 
required  for  subsequent  growth  of  an  amount  a  On 
this  b  asis,  it  can  be  expected  in  some  problems  that 
t,  will  be  negligible  compared  with  the  total  time 
required  for  failure  of  a  structure 

Whether  or  not  T  equals  T,  depends  on  the  he- 
reditarv  and  multiaxial  stress  characteristics  of  the 
failure  rone,  the  physical  environment  (since  the  left 
edge.  ;  =  o.  is  at  the  surface  during  0  <  t  <  (,)  and 


on  the  state  of  the  initial  crack  tip  In  an\  event  i! 
should  not  be  assumed  that  these  energies  ate  equal 
without  adequate  support  from  experimental  data 

4  3  Energy  Release  Rate 

The  7.  integral  is  equal  to  the  energv  release  rate 
for  self-similar  crack  growth  (An*  =  An?i2l  over 
Os  |  s  q)  in  the  reference  elastic  problem,  w  hich  is  a 
well-known  and  useful  relationship  for  experimental 
and  theoretical  determination  of  7,  Namely. 

J,=  -dP„'dA  (311 

where  dP,  is  the  change  in  potential  energy  of  the 
reference  elastic  continuum  and  applied  forces.  d.A 
is  the  increase  in  crack-plane  area  This  relationship 
is  easily  derived  by  considering  the  work  done  hv  a 
general  three-dimensional  elastic  continuum  on  a 
failure  zone  which  is  advanced  an  amount  da  along 
a  segment  of  the  crack  edge  (Schapery  1 9*4) 


5.  Power  Law  Nonlinearity 

Here,  the  specific  effects  of  nonlinearity  on  the  open¬ 
ing  mode  of  crack  growth  are  illustrated  using  a 
power  law  nonlinear  material  By  definition,  the 
pseudo  strain  energv  density  in  this  case  for  isotropic 
or  anisotropic  materials  is  a  homogeneous  function 
of  degree  .V  +  1 . 

<t>(ccf)  =  :,c')'*‘<J>(rf )  |32) 

where  ,Y  and  c  are  constants  and  1  denotes  absolute 
value.  For  notanonal  simplicity,  the  superscripts  <• 
and  R  will  be  omitted  in  succeeding  equation'  until 
specific  results  are  applied  to  viscoelastic  crack 
growth  analysis  With  the  definitions 

f,'=rr,.  <t>'  =  <t>(r  ).  sgn(r)  =  sign  of  c  (33) 

it  follows  from  Eqns  (15)  and  (32)  and  application 
of  the  chain  rule  to  h<t>  He,,  that 

c<<l>  ft  9> 

—  =  sgn(f)lc[  : — .  o,  =  sentr  ).c!  V  (34i 

iff,  rtf, 

For  a  uniaxial  strain  state  the  second  relation  implies 
<7n  —  sgn(f,i)ifu!';  the  same  form  is  predicted  for 
uniaxial  stress  (Schapery  1981)  Equation  (32l  con¬ 
tains  as  a  special  case  the  power  law  strain  energv 
function  commonly  used  for  isotropic  materials  (  Rice 
and  Rosengren  19681 

Let  us  introduce  a  set  of  dimensionless  coordinates 
and  mechanical  state  variables. 


=  sgni  o,. )  o„  om' s  r„ 

it  *  sgn(om).o,  om  1  'u  n  (35i 

The  quantities  and  a  are  parameters  with  dimen¬ 
sions  of  stress  They  are  independent  of  x  .  and  are 
introduced  for  the  purpose  of  defining  a  dimen- 
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sionless  failure  stress  distribution  / and  strain  energy 
density  <1\. 


f  =  <>,  rrm.  <tv  ®  <t>  o„  (36) 


The  parameter  n„  is  assumed  to  depend  at  most  on 
aging  time  When  these  definitions  are  substituted 
into  Eqns  ( 10).  (11)  and  ( 15).  and  Eqn.  (34)  is  used, 
the  field  equations  become 


S&„ 

St 


t. 


'( 


dii,  Su 
(it,  +  St, 


~st~ 


(37) 


w  here  <t\  ^  <t>\(  f,.l  Allowance  is  made  for  negative 
valuer  of  it,  through  the  factor  om.  for  unloading  of 
a  viscoelastic  material  this  stress  could  be  neeative 
•n  limited  periods  of  time  vs  ithout  crack  face  contact 
Locating  the  origin  of  the  coordinates  at  the  crack 
tip  P  (n  =  tl  tn  Fig  2).  and  assuming  the  crack  faces 
are  traction  free  outside  of  the  failure  zone.  then, 
with  the  additional  definition  q  5=5  (2  n)  m  -a,,  the 
solutions  to  Eqn  (37)  must  satisfv  the  traction  bound¬ 
ary  condition  d;;  =  /(»))W(  I  -  rj)  for  r;  >  0  as  well  as 
meet  the  symmet rv  conditions  associated  with  the 
opening  mode  These  solutions  will  depend  at  most 
on  the  dimensionless  coordinates  i,.  apart  from  pa¬ 
rameters  in  <l\  and  f  and  those  arising  from  con¬ 
ditions  imposed  on  the  continuum  surrounding  the 
failure  zone  For  example,  the  solution  for  dimen¬ 
sionless  crack-opening  displacement  is  simplv  Aii;  = 
gt  r;i  f  rom  Eqn  (35). 


=  sgn(  evict  ev  oT  '  g(  i/ 1 

(38) 

Equations 

( 18).  (20).  (35)  and  (38)  yield 

1 1’- : ! '  A 

a  -  j  —  ,  - 

1  tV  \Om'l  ( 

( 39a ) 

where 

( ’  rfg 

A  =  )  /('/'^dt; 

(39b) 

From  Eqns  (38)  and  (39). 

A  It;  =  (A.  <v/,)g(  7/ ) 

(40) 

The  results  in  Eqns  (38-40)  are  not  limited  to  a 
small-scale  failure  zone  and  are  valid  for  materials 
which  age  and  are  physically  nonhomogeneous  with 
respect  to  z;;  the  dependence  of  <t>  on  x:  and  t  and 
related  material  parameters  has  been  suppressed  to 
simplify  the  notation  The  dimensionless  displace¬ 
ment  g.  stress / and  integral  I,  could  in  general  depend 
on  all  dimensionless  parameters  that  influence  con¬ 
tinuum  deformation  in  the  neighborhood  of  the  crack 
tip  Considerable  simplification  in  this  dependency 
results  if  the  continuum  is  physically  homogeneous 
(at  least  locally)  and  if  a  is  small  relative  to  the 
distance  to  other  geometrical  features  The  immedi¬ 
ate  continuum  surrounding  the  crack  tip  neigh¬ 
borhood  can  then  be  analyzed  by  using  a  method 


similar  to  that  in  Rice  and  Rosengren  (196.8)  for  an 
isotropic  power  law  material  With  a  polai  coordinate 
system  (r.  6 )  centered  at  P  (compare  with  Fig  2)  it 
is  found  that  for  r  $>  q, 

u,  =  r ( /, '/)'  ,wl7i,(») 

t„  =U,/r)'  "•'/(«) 

=  (/./f)',N*  "g,A0)  (41) 

wherer=r  n  The  integral  /,  is  related  to  A  throueh 
Eqn  (39a).  and  the  functions  of  6  are  determined 
from  the  solution  of  Eqn  (37)  and  the  condition  of 
traction-free  crack  surfaces  for  f  >  0.  The  failure 
zone  traction  on  the  crack  plane  is  neglected  in  solv¬ 
ing  this  boundarv  value  problem;  its  effect  enters  the 
solution  in  the  form  of  /,  after  path  independence  of 
J ,  is  taken  into  account  through  Eqn.  (39a)  (which  is 
based  on  Eqns  (18)  and  ( 2U ) )  The  dependence  of 
Eqn  (41 )  solutions  on  r  is  the  same  as  found  hv  Rice 
and  Rosengren  ‘  1968).  but  the  l>  dependence  is  not 
necessarily  the  same  because  the  strain  energv 
density  used  here  is  more  general  Such  detailed 
information  is  not  needed  in  this  studv.  ralhet.  it  is 
sufficient  to  know  that  the  entire  influence  of  external 
loads  and  other  far-field  parameters  and  geometry 
(such  as  other  cracks)  on  the  dimensionless  mech¬ 
anical  state  variables  in  Eqn  (35)  in  the  neigh¬ 
borhood  of  the  crack  tip  is  accounted  tor  hv  /. 

With  the  small-scale  failure  zone,  given  the  shape- 
function  fit))  for  the  failure  stress  distribution  and 
not  using  Eqn  (39b).  it  may  he  concluded  that  the 
predicted  shape  function  g(i)l  for  the  crack-opening 
displacement  can  depend  at  most  on  /,.  ,V  and  aging 
time  r  Equation  (39b)  is  thus  an  implicit  relationship 
for  predicting  /,.  The  significance  of  Eqn  (34h!  may 
be  clarified  by  imagining  an  analysis  in  which  external 
loading  is  increased  from  zero  by  using  increasing 
values  of  /,  in  Eqn  (41).  but  not  using  Eqn  (39bl. 
until  infinite  stresses  and  strains  at  P  arc  removed 
For  the  linear  problem.  ,N  =  I .  it  can  be  verified  that 
the  removal  occurs  when  /,  is  equal  to  the  integral  in 
Eqn.  (39b);  with  this  value  for  /,.  Eqn  (39a)  becomes 
identical  to  a  dimensionless  form  of  the  Uurenbiatt 
(1962)  condition  for  finite  crack  tip  stress  after  the 
familiar  relationship 

A  =  (1  -  v:)K\  £«  (42) 

is  employed  (where  A',  is  the  opening-mode  stress- 
intensity  factor.  £r  is  Young's  modulus  and  i-  is 
Poisson's  ratio)  Recall  that  Eqn  (39)  comes  from 
Eqns  (18)  and  (19)  (through  the  special  case  of  Eqn 
(20))  in  which  the  portion  of  the  contour  C-  around 
the  tip  P  vanishes  on  account  of  the  physical  require¬ 
ment  of  finite  stresses  and  strains  Thus,  it  appears 
that  Eqn  ( 18).  in  which  Eqn  (20)  is  used  for  J,.  can 
be  interpreted  as  a  generalization  of  Barenblatt  s 
condition  to  nonlinear,  anisotropic,  viscoelastic 
media  with  an  arbitrarily  large  failure  zone  For 
such  a  general  situation,  these  equations  mnv  not 
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guarantee  that  the  crack  faces  do  not  overlap,  and 
thus  additional  analysis  involving  contact  phenomena 
could  he  needed 

It  hav  been  shown  that  /,  is  a  constant  when  /  = 
/(i/l.  except  for  possible  aging  It  should  be  clear 
that  this  is  true  even  when  /  depends  also  on  the 
displacement  shape  function  g  In  this  latter  case  one 
must  solve  simultaneous!)  for  /.  g  and  /,  using  the 
dimensionless  equations. 

5  I  Cnuk  Growth 

The  analysis  in  this  section  has  so  far  been  limited  to 
a  studs  of  the  reference  elastic  problem  With  the 
addition  of  the  superscript  R  to  strains  and  displace¬ 
ments.  thev  mav  be  used  with  the  results  in  Sect  4 
Equations  ( 23 )  and  (24)  for  fracture  initiation  were 
derived  assuming  o,  does  not  vary  with  represented 
here  bv  t  =  1  since  o„  =  o,  In  this  case  the  equations 
for  initiation  time  t,  are  very  simple  and  are  not 
limited  to  a  small-scale  failure  zone  If  df  d>i  *  0  it 
is  helpful  to  use  the  results  for  a  power  law  material 
to  predict  t  from  Eqns  (40)  and  ( 14)  we  obtain  the 
pseudo  and  viscoelastic  displacements  (with  g{  I )  = 
g  t. 

An ",  =  J,g,  ,ij,.  Au:,  =  {DdAu ?„}  (43) 

1  1  "  - .  ■  ■  ■  -  v  .hi  Of  defined  as  o.  at  2  -  n  w  ithout 

-  •  •  •’ i* v  the  analvsis  If  aging  is  negligible, 
t”.  '  ■  s  .  ■  ■ . ;  i  -p  m  f  qn  I  22  i  reduces  to  Eqn  (23 1 
e»  •  I  w  •'•o.Jified  bv  the  constant  factor  I,  g , 

f  1 '  '  2,J  1  li’t  i  rack  speed  depends  on  the  si/e 

1  •  ’  .••.  ’  >rte  I  his  si/e  now  assumed  small 

1  '  •  •  •'<•*  i  i w  material  is  given  bv  Eqn  (3yji 

1  -  's  ■  si .  il  sicnilis.m.e  of  this  result  First . 

"  ■  :  "•-•••  ’  mav  be  interpreted  as  a  v icld ‘tress 

'  •  raaiii  when  \  is  small  using  <i>  =  a  <J\ 

'  •  wh  I  qn  (32)  and  one-dtmenstonal  or 

P'  p  v:  't’.ii  loading  thri'ugh  a  stress  then  r'  - 
•  ■  heme  the  strain  is  small  when  o  <  o,  and 

.  w  bt  ••  o  (t  d  \  .*  1  Also.  i».  mav  be  con- 
si  'lie  m'rinsie  strength  if  it  is  taken  to  be  the 

r  o  m..m  v.,iue  (with  respect  to  location  2i  of  the 
tii  ,'e  Z  'ne  stress  durinc  crack  growth  Thus,  the 
si.'-,  a  m  Eqn  (  Vuj  |  increases  with  increasing  yield 
s"e-'  tits!  is  sensitive  to  the  rativ’  of  vield  stress  to 
in'tin-i.  sttencth.  given  J  In  elasioplastic  fracture 
mt'sft. lines  n  is  common  practice  to  assume  ar  is 
proportional  to  the  v  leld  stress  and  f  =  1 .  in  this  case, 
the  standard  result  a  ~  J ,  o.  is  recovered  since  l<  is 
a  const,  i  ni  More  general  I  v .  o_  mav  varv  with  n  (and 
o’her  p.irame(ers)  so  that  a  J .  is  not  necessarily 

C'  ' I' -I  111! 

I  he  results  so  far  have  been  expressed  in  terms  of 
material  functions  for  the  failure  zone,  a ,  and  T 
I  hese  quantities  may  depend  on  a  and  other  variables 
suMi  as  temperature  and  moisture,  especially  if  the 
zone  is  viscoelastic  One  could  explicitly  incorporate 
specific  phvsical  models  of  the  failure  zone  in  order 
to  complete  the  formulation  of  the  theory  Assuming 


the  failure  process  for  material  elements  on  the  crack 
plane  is  unaffected  bv  stresses  prior  to  the  arroa1  or 
near  arrival  of  the  crack  tip.  such  an  analvsis  lor  a 
small-scale  failure  zone  would  yield  a  function  «  = 
a{JJ.  with  possible  dependence  on  phvsic.il  par 
ameters  and  material  constants  such  as  icnipcratuic 
age  and  nonlinear  exponent  \  .  that  the  instant, me- no 
2,  (hut  not  past  values)  determines  the  instantaneous 
a  follows  from  the  fact  that  the  crack  tip  neigh¬ 
borhood  is  surrounded  by  a  stress  field  which  is 
defined  solely  by  the  current  value  of  J.  (i  e  .  Eqn 
(41)  but  in  terms  of  the  dimensional  phvsical 
variables),  assuming  7,  changes  slowly  enough  that  a 
is  essentiallv  constant  during  the  failure  time,  n  a 
Thus,  when  the  dependence  of  a  on  basic  material 
parameters  is  not  of  direct  interest,  one  could  charac¬ 
terize  fracture  behavior  with  the  function  a  -  nU  ) 
found  directly  from  experimental  data  II  theoretical 
investications  indicate  a  range  exists  for  which  ./<. 
dl  <  0  unstable  crack  growth  (  possibly  in  steps  of 
stop-start  behavior!  may  result  and  the  intrinsic 
a  -  J,  function  for  slow,  continuous  gross t h  could 
not  be  found  experimentallv  Indeed,  for  this  case 
an  appropriate  analysis  may  involve  initiation  and 
dynamic  arrest  phenomena 

5  2  Power  Law  Viscoelasticu\ 

The  creep  compliance  for  viscoelastic  materials  can 
often  be  expressed  as  a  power  law  m  time  Ihus.  as 
an  important  special  case,  consider  nonaging  material 
behavior  for  which  D(t  -  r)  =  Ddt  -  ri"  t where  D 
and  n  are  positive  constants)  in  order  to  predict  the 
form  of  the  equation  a  =  a(J,)  using  certain  specified 
failure  zone  characteristics  A  small-scale  failure 
zone  with  a  shape  factor  f  e  a,  o„  that  depends  at 
most  on  g  and  i/  =  J  a  is  assumed  In  general  r>„ 
and  T  may  depend  on  crack  speed  and  other  par¬ 
ameters  If.  however.  ct„  and  E  are  constant  (so  that 
this  zone  exhibits  elastic-like  behavior)  Eqns  (24) 
and  (39ai  yield 

a  -  J\  (44 | 

where  p  =  1  +  ( 1  n )  If  a  (instead  of  o m)  and  V  are 
constant.  Eqn  (29)  predicts  p  =  I  n  The  failure 
zone  dimensions  are  defined  by  a  and  crack-opening 
displacement  Au-  at  e  =  a.  if  these  two  quantities 
are  constant,  one  may  use  Eqn  (39a)  to  solve 
for  om.  and  then  use  Eqn  (2?)  for  Att:  and  Eqn 
(40)  for  Au?  at  17  =  *  a  =  1  to  find  that  p  = 
[n(l  •+■  ,V)j  1  Notice  that  only  in  this  last  case  does 
the  nonlinearity  exponent  appear  in  the  equation  for 
P 

Experimentally  determined  values  of p.  when  com¬ 
pared  with  the  results  for  these  and  other  cases,  mav 
be  helpful  in  determining  failure  zone  characteristics 
and  guiding  theoretical  model  development  The 
values  of  p  =  3  and  n  =  0  5  were  obtained  from 
experimental  data  on  a  g I ob alls  linear  elastomer  this 
is  consistent  with  the  assumption  of  constant  E  and 
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o.,  ( Schapery  19~5pp  549-62)  Some  studies  of  crack 
growth  m  mctaK  undergoing  secondare  creep  li  e  . 
a  =  I  )  pi. mde  values  of  p  =  1  (Landes  and  Beglev 
19'(i.  Nikbin  ct  al  1976).  since  A  <t  I  tvpicnllv .  the 
data  are  consistent  » it h  the  last  two  cases  for  constant 
values  ot  I  and  o  or  An-,,  and  a  It  should  he  added 
that  the  exponent  .V  used  here  pertains  to  nonlinear 
behavior  of  the  continuum  in  the  neighborhood  of 
the  crack  tip  The  relationships  are  valid  even  if 
nonlinearity  of  material  far  from  the  tip  (relative  to 
a )  does  not  obey  a  power  law  .  or  is  linear  (i.e. ,  A  = 
II  w  hile  its  local  behav  lor  is  described  using  X  *  1. 
The  remote  nonlinear  behavior  does  not  affect  the 
e  -  7,  equation,  but  it  does  affect  the  functional 
dependence  of 7,  on  the  externally  applied  loads 

6  Com  Installs 

Through  the  use  of  a  nonlinear  single-integral  consti- 
tutive  equation,  generalized  ./-integral  theory  and 
certain  approximations,  relatively  simple  crack 
growth  equations  based  on  mechanical  work  for  time 
of  initiation.  Eqns  (23)  and  (24).  and  speed.  Eqn. 
(29)  were  derived  Onlv  the  opening  mode  of  growth 
has  been  analyzed  in  detail  here,  but  the  basic 
relationships  m  Sects  1-3  could  be  used  for  shearing 
and  mixed-mode  conditions  Rheological  properties 
of  the  continuum  are  reflected  in  the  creep  com¬ 
pliance  D  and  the  generalized  integral  7.  The  length 
a  of  the  zone  of  failing  material  at  the  crack  tip 
appears  in  the  equation  for  speed,  and  it  provides  a 
scale  which  determines,  in  effect,  the  magnitude  of 
local  strain  rates  resulting  from  crack  growth  In 
general  n  is  not  constant,  instead,  it  is  related  to  7. 
through  Eqns  (IS)  and  (20).  which  reduce  to  Eqn 
|5-<I  with  power  law  nonlinearity 

All  of  these  results  come  from  an  analysis  of  the 
continuum  in  the  vicinity  of  the  crack  tip  Thus,  for 
limited  crack  growth  it  is  not  necessary  for  the  entire 
body  to  he  represented  by  Eqn  (1)  and  the  local 
creep  compliance  A)  However,  if  it  is.  Eqn  (31)  may 
be  employed  to  predict  7,  with  an  arbitrary  amount 
of  crack  grow  th  bv  using  theoretical  and  experimental 
methods  similar  to  those  already  developed  for  time- 
independent  materials  In  this  expression  E\  and  7, 
are  defined  just  as  for  elastic  materials,  but  dis¬ 
placements  and  strains  are  pseudo  quantities,  uR  = 
IE  7u,)and  rR  =  r'  =  (£  de„],  respectively.  Stresses 
and  loads  are  the  actual  physical  quantities  For 
example,  experimentally  measured  displacements  u, 
would  be  converted  to  pseudo  displacements  u,R 
through  Eqn  (K|  before  evaluating  7,  from  lest  data. 
The  coefficient  £„  is  a  free  constant,  and  could  be 
taken  as  the  modulus  or  reciprocal  of  compliance  at 
some  specified  time 

The  bibhographv  includes  publications  which 
devenbe  various  models  for  characterizing  and  pre¬ 
dicting  crack  growth  in  different  types  of  materials 
Sev  eral  models  are  special  cases  of  the  general  theory 


presented  here  For  secondary  creep  li  e  .  nonlinear 
viscous  behavior),  experimental  data  on  crack  speed 
have  been  successfully  correlated  in  terms  of  a  C  ‘ 
parameter  by  many  investigators  (e  g  .  Landes  and 
Begley  1976);  the  function  a  -  a(C’)  typicallv  obev s 
a  power  law  This  characterization,  whether  or  not 
the  power  law  is  obeyed,  is  obtained  from  the  present 
theory  by  using  the  relaxation  modulus  E  = 
£R6(r-  r)  in  Eqn.  (8).  which  reduces  the  pseudo 
variables  to  strain  rates  (f,r  =  c ,,)  and  velocities 
(u,R  =  u,);  in  turn,  7.  becomes  C*  The  use  here  of 
nonsingular  strains  at  the  crack  tip,  in  contrast  to 
earlier  work  on  secondary  creep,  leads  to  a  simple 
physical  interpretation  of  7.  for  viscous  media  Spe¬ 
cifically.  from  Eqn.  (29)  with  k  =  1/3  and  !)  = 
£r  1  ( z  —  r) .  corresponding  to  the  aforestated 
modulus,  we  find  7.  =  3(2Ta  a)  =  3  x  (mechanical 
power  input  to  the  failure  zone  per  unit  areal  f  or 
viscoelastic  materials  7,  has  a  simple  physical  mean¬ 
ing  if  the  continuum  is  essentially  elastic  except  for 
a  small  amount  of  material  around  the  crack  tip  In 
this  case  7,  is  approximated  the  energy  release  rate. 
-dP/dA,  since  u,R  -  u;.  aiso.  T  £„£>  in  Eqn  (29) 
may  be  properly  called  an  effective  fracture  energv 
Crack  growth  in  homogeneous  linear  viscoelastic 
media  has  been  characterized  traditionally  in  terms 
of  stress-intensity  factors  because  of  the  many  situa¬ 
tions  in  which  viscoelastic  stresses  are  the  same  as  in 
elastic  materials  With  isotropy .  linearity .  a  constant 
Poisson's  ratio,  homogeneity  and  locally  plane  strain , 
it  is  found  that  Eqn  (41 )  (after  expressing  it  in  terms 
of  the  stress-intensity  factor)  and  Eqn  (39a  I  yield 
Eqn  (42)  Substitution  of  Eqn  (42)  into  Eqns  (24) 
and  (29)  results  in  the  same  crack  growth  relations 
derived  earlier  (Schapery  1975).  Similar  results  in 
terms  of  stress-intensity  factors  for  linear  behavior 
exist  without  the  restrictions  of  constant  Poisson  s 
ratio,  isotropy  and  homogeneity  (Brockway  and 
Schapery  1978,  Schapery  1975.  1 978 1  This  obser¬ 
vation  includes  cracks  between  dissimilar  media  (i  e  . 
adhesive  fracture)  if  both  materials  are  incom¬ 
pressible  or  one  is  relatively  rigid  and  the  other  is 
incompressible;  without  the  incompressible 
behavior,  tensile  and  shearing  stresses  act  simul¬ 
taneously  on  the  crack  plane,  and  a  complex  mixed- 
mode  condition  generally  exists  at  the  crack  tip  The 
nonlinear  theory  in  Sect  4  is  applicable  to  adhesive 
fracture  if  the  mixed-mode  state  does  not  exist  and 
D  is  the  same  for  both  materials  (or  one  is  relatively 
rigid);  T,  and  T  in  Eqns.  (24)  and  (29)  obviously  have 
to  be  interpreted  as  the  fracture  energies  for  the 
particular  material  combinations  involved 

Primary  creep  of  metals  and  ceramics  is  cus¬ 
tomarily  represented  by  using  a  nonlinear  power  law 
model  of  viscous  flow  with  strain  hardening  When 
crack  growth  exists  this  behavior  is  not  described  bv 
Eqn  (1).  and  the  correspondence  principle.  Eqn 
(13).  does  not  apply  However,  with  proportion.il 
loading  and  a  stationary  crack,  the  strain-hardening 
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constitutive  equation  becomes  identical  to  that  of  a 
nonlinear,  aping  elastic  body  (e  g  .  Riedel  1981) 
liquation  ( I )  takes  the  appropriate  form  by  using 
I  Kl)  =  /i(M  A  suitable  choice  of  /i(i)  and  a  power 
law  potential  <J*(r ' ).  where  the  only  time  dependence 
is  through  rf.  =  ejh.  yields  the  desired  primary 
creep  behavior  as  well  as  a  singular  stress  field 
parameter  ('!  (Riedel  1981)  which  is  identical  tod,. 
For  a  small-scale  failure  zone,  the  early  stages  of 
crack  growth  are  determined  by  the  stress  field  sur¬ 
rounding  a  stationary  crack  tip  and  therefore  by  7, 
However,  a  theory  does  not  appear  to  be  available 
for  relating  crack  growth  to  a  single  parameter  when 
the  amount  of  growth  becomes  comparable  to  the 
scale  of  the  initial  singular  stress  field 

It  should  be  observed  that  a  different  type  of 
primary  creep,  such  as  that  used  to  characterize 
pol\  mers.  is  contained  in  Eqn  ( 1 )  through  the  depen¬ 
dence  of  creep  compliance  /)  on  t  -  r.  Also,  through 
the  use  of  aging  time  in  <t>  =  4>(f ' .  r).  one  may 
characterize  different  types  of  nonlinear  behavior  at 
short  and  long  times  even  though  the  creep  com¬ 
pliance  itself  is  not  a  function  of  stress  Thus,  in  spite 
of  the  apparent  simplicity  of  Eqn.  ( 1 )  compared  with 
other  available  nonlinear  viscoelastic  constitutive 
equations,  it  is  really  quite  general  and  yet  permits 
the  use  of  the  correspondence  principle  with  large 
amounts  of  crack  growth  Nevertheless,  the  accuracy 
of  Eqn  ( I )  for  many  materials  under  various  stress 
histories  is  not  yet  established 

The  deformation  and  fracture  analysis  has  been 
formulated  assuming  small  strains  and  rotations  for 
the  continuum  This  restriction  is  needed  because  the 
correspondence  principle  is  not  generallv  valid  with 
nonlinear  strain-displacement  relations  An  approxi¬ 
mate  theory  with  large  deformations  is  described  bv 
Schapery  (19891  Another  possible  complication  not 
explicitly  treated  here  is  crack  tip  heating,  which 
can  be  very  significant  in  polymers  because  of  their 
tvpicallv  low  thermal  conductivity  and  high  ultimate 
strains  If  this  heating  is  sufficiently  localized,  it  could 
be  included  in  the  characteristics  of  the  failure  zone, 
and  thus  not  complicate  the  continuum  analvsis 

T  he  idealization  of  a  slender  failure  zone  is  quite 
realistic  for  crack  tips  in  many  materials,  as  well  as 
for  craze  tips  in  plastics  In  considering  the  zone 
shape  and  size,  it  should  be  recalled  that  only  the 
part  of  the  body  which  does  not  obey  Eqn.  (15)  has 
to  be  included  in  the  failure  zone.  In  some  cases  it  is 
helpful  to  represent  part  of  the  failing  material  (such 
as  craze  I  as  a  surface  loading,  consequently  account¬ 
ing  for  it  in  7,  As  a  further  generalization  it  is  shown 
bv  Schaperv  (1984)  that  important  aspects  of  the 
theorv  remain  valid  even  if  distributed  damage,  such 
as  microcracks,  is  included  in  the  mathematical 
model  of  the  continuum 

Considerable  progress  has  been  made  in  the  last 
twenty  vears  in  understanding  effects  of  time-depen¬ 
dent  rheological  properties  on  crack  growth.  How- 
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ever,  most  basic  experimental  and  theoretical  inves¬ 
tigations  have  been  for  the  opening  mode  1 4  cohesiv e 
fracture  in  essenliallv  homogeneous  materials  In  this 
same  time  period,  use  of  liber-reinforced  materials, 
especiallv  plastics,  in  load-bearing  structures  has 
increased  dramaticallv  Thus,  considering  the  com¬ 
plex  interactions  which  may  occur  at  various  scales 
and  the  importance  of  interfaces  in  fibrous  lamina 
and  laminates  and  in  other  multiphase  materials 
there  is  an  especially  important  need  for  more  basic- 
studies  on  mixed -mode  crack  growth  adhesive  crack¬ 
ing  and  the  effect  of  nonhomogeneous  time-depen¬ 
dent  properties  on  all  types  of  crack  growth 

See  also  Time-Dependent  Fracture  Mechanisms.  Mech¬ 
anics  of  Materials:  An  Overview 
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stresses  that  arc  onlv  a  small  lijstion  of  the  -  'hoi 
strength,  this  condition  is  rttet  at  vu.li  large  err.  • 
levels  of  the  cracked  svsicrtt  that  thermal  mop 
cannot  play  a  significant  role  in  sot'ci”"L  r 
response  Hence  the  onlv  rate-dependent  ft.r.  •.;• 
in  elastic  svstems  at  low  temperatures  arc  "i  ' 
stress-corrosion  cracking  or  statu  !af;guc  tv  pc  »(>•. 
slow  crack  growth  under  subcriiK.il  von.!:'i 
becomes  possible  bv  rate  processes  .it  via.k  tip  - 
rosion  or  contamination  These  '  v  pcs  of  j.iik  gr  ow 
are  discussed  elsewhere  tsee  %."«••»  <  ■•rrou.ot  ( 
trig  Suh,  rut,  al  (  nn  k  GVow/'i  in  <  <  ’,imn 

The  time-dependent  fractures  t  modt  'ed  fieri1 
those  which  involve  time-dependent  qu.i-i  fr 
eeneous  damage  jwuiiu'Ci  u:  ..n  let  s'tr 
microscopic alls  mhom.'gene  'inlv  sro-.  .  J  u’  ‘ 
homogcnc'vislv  defotmuig  v  I  ••  p 
not  son; .nit  initial  .  r  a.  ►  -  ■  "  '  11  -  ■  ■  ■ 

sulvnii.  .,1  si/c  i'h  rv...  •••.-.  i  ■  ■ 
accumulation  (>t  damage-  -tits'  :  .e  lot’-  " 

more  concentrated  atounJ  sites  •>(  p-c’-  ,  ■.;■ 
or  at  tips  ol  subs  ritual  sinks  |  p(-  tin.n  .ml  a 
usuallv  tune  dependent  growtft  ot  the  i.x1 
takes  comparative!'  little  time  In  most  ot 
processev  the  time  dependeiist  largclv  reside  -  u  ' 
deformation  processes  whi.tr  vbronologi  .R . 
the  mivroprocesses  that  huil.l  up  t tie  mis  •  ■  . 

ccntrations  even  in  the  pttvcn.e  of  o'lict  snui 
processev  that  tr\  to  dissipate  these  same  stress  , 
centrations  Onlv  rarely  does  the  time  dependtr 
reside  in  true  processes  of  decohc'ion  on  the  m 
ecular  svjle  The  phases  of  the  time-dependent  It 
lure  that  involves  crack  grow  th  in  a  creeping  solid 
dealt  with  elsewhere  (see  1  one  Dependent  f>a,t;, 
Continuum  Aspects  ol  Crock  drouth I 


/  Forms  of  Damage 
I  I  Rupture  versus  Damage 

A  bat  in  tension,  as  shown  in  Fig.  1.  capable  ot 
undergoing  time-dependent  deformation  in  steads 
state  under  a  constant  stress  at  h'gh  tempera! urc  will 
rupture  as  a  result  of  gradually  accelerating  deform¬ 
ation  when  under  a  constant  applied  load  as  the  area 
continuously  decreases  If  the  steady  state  strain  rate 
f  is  a  power  function  of  tensile  stress  o.  in  the  form 
f  =  fl(T)(  o/On)”'.  where  B{  T).  cv,  and  m  are  material 
parameters,  the  time  to  complete  rupture  is 


Time-Dependent  Fracture:  Mechanisms 

Fracture  of  solids  is  accomplished  by  the  formation 
and  growth  of  cracks  under  stress  In  purely  elastic 
bodies  the  impending  propagation  of  a  crack  is  associ¬ 
ated  with  a  "sxstem"  instability  in  which  ever  larger 
amounts  of  elastic  strain  energy  are  released  by  a 
growing  crack  than  is  required  to  produce  the  sur¬ 
faces  of  the  cracks  Under  usual  conditions  of  applied 


f,  =  1/mf,  ( I  ) 

where  f,(  =  fl(  7  Ko,  an)-")  is  the  initial  creep  strain 
rate  As  Fig  2  shows  for  two  typical  cases  of  m  =  5 
and  m  =  8  for  engineering  alloys,  the  acceleration  of 
creep  in  the  final  stages  is  extremelv  rapid 

In  the  maiority  of  cases  where  steady-state  deform¬ 
ation  is  attainable  only  after  a  pnmarv  creep  transi¬ 
ent,  the  overall  extension  of  the  bar  is  somewhat 
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ABSTRACT 

In  most  materials,  macrocrack  extension  is  accompanied  by  inelastic 
phenomena  (such  as  microcracking  or  plastic  deformation)  throughout  a  region 
surrounding  the  crack  tip.  Immediately  ahead  of  the  crack  tip,  strain 
localization  occurs  in  a  small  volume  of  heavily  damaged  material  referred  to 
as  the  failure  zone  or  fracture  process  zone.  In  this  study,  the  failure  zone 
and  the  surrounding  zone  of  inelastic  material  are  treated  as  two  distinct 

regions.  The  failure  zone  is  assumed  to  be  thin  relative  to  its  length  and  is 
represented  in  a  two-dimensional  finite  element  model  as  tractions  which  act 
across  the  crack  faces  near  the  tip.  An  opening  mode  of  crack  tip  deformation 
is  assumed.  The  normal  traction  at  any  point  on  the  crack  surface  in  the 
failure  zone  is  specified  as  a  decreasing  function  of  the  crack  opening 
displacement  which  vanishes  after  a  critical  value  of  displacement  is 

reached.  Two  different  rate-independent,  inelastic  continuum  characteri¬ 
zations  are  used;  one  models  metal  plasticity  and  another  represents 

microcracking  in  brittle  materials.  Both  constitutive  models  allow  for  the 
definition  of  a  generalized  J-integral  developed  by  Schapery,  which  has  the 

same  value  for  most  paths  around  the  crack  tip  for  realistic  distributions  of 
plasticity  or  damage  in  the  material  surrounding  a  stationary  or  propagating 


4 

4 

crack.  Tnis  path  independence  is  verified  numerical/  for  a  c^ack  growing 
under  conditions  of  small-scale  inelasticity,  and  the  equivalence  between  J 
and  the  work  input  to  the  last  ligament  of  material  in  the  failure  zone  is 
demonstrated.  Steady-state  crack  growth  is  stuided  in  two  different  specimen 
geometries.  Simplified  J  integral  analyses  are  used  to  estimate  the  work 


input  to  the  failure  zone  for  these  steady-state  problems.  The  J  integral 
estimations  are  compared  with  finite  element  results  to  determine  the  accuracy 
of  the  simplified  analyses. 
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Techniques  for  Characterizing  Damage  Zones  in  Composite  Materials 

by 

Richard  Dale  Tonda 

Interdisciplinary  Engineering,  August  1987 
ABSTRACT 

Experimental  methods  for  determining  the  existence  of  a  work  potential, 
useful  in  characterizing  material  response  even  in  the  presence  of  significant 
damage,  are  developed  and  described.  The  underlying  fracture  theory  and 
motivation  for  the  development  is  briefly  discussed.  The  theory  involving  the 
work  potential  and  the  existence  and  use  of  a  damage  parameter  is  reviewed,  as 
is  the  basic  composite  laminate  theory  and  underlying  constitutive 
formulation.  An  experimental  program  which  addresses  these  factors  through 
the  biaxial  loading  of  composite  tubes  is  presented.  The  results  are 
analyzed,  incorporating  viscoelasticity  theory  for  power  law  time  dependence 
of  the  matrix  material,  and  it  is  demonstrated  that  the  theory  provides  a  less 
complicated  presentation  of  the  data  by  eliminating  the  superimposed  effects 
of  time  dependent  behavior. 


PH.D.  DISSERTATION 


Finite  Element  Analysis  of  Crack 
Growth  in  Inelastic  Media 
by 

Joe  Randall  Weatherby 
Mechanical  Engineering,  May  1986. 

Co-Chairmen  of  Advisory  Committee:  Dr.  R.A.  Schapery 

Dr.  W.L.  Bradley 

ABSTRACT 

In  most  materials,  macrocrack  extension  is  accompanied  by  inelastic 
phenomena  (such  as  microcracking  or  plastic  deformation)  throughout  a  region 
surrounding  the  crack  tip.  Immediately  ahead  of  the  crack  tip,  strain 
localization  occurs  in  a  small  volume  of  heavily  damaged  material  referred  to 
as  the  failure  zone  or  fracture  process  zone.  In  this  study,  the  failure  zone 
and  the  surrounding  zone  of  inelastic  material  are  treated  as  two  distinct 
regions.  The  failure  zone  is  assumed  to  be  thin  relative  to  its  length  and  is 
represented  in  a  two-dimensional  finite  element  model  as  tractions  which  act 
across  the  crack  faces  near  the  tip.  An  opening  mode  of  crack  tip  deformation 
is  assumed.  The  normal  traction  at  any  point  on  the  crack  surface  in  the 
failure  zone  is  specified  as  a  decreasing  function  of  the  crack  opening 
displacement  which  vanishes  after  a  critical  value  of  displacement  is 
reached.  Two  different  rate-independent,  inelastic  continuum  characteri¬ 
zations  are  used;  one  models  metal  plasticity  and  another  represents 
microcracking  in  brittle  materials.  Both  constitutive  models  allow  for  the 
definition  of  a  generalized  J-integral  developed  by  Schapery,  which  has  the 
same  value  for  most  paths  around  the  crack  tip  for  realistic  distributions  of 
plasticity  or  damage  in  the  material  surrounding  a  stationary  or  propagating 


vr 
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